LOCALIZED WAVES: A NOT-SO-SHORT REVIEW (t) 



Erasmo Recami 

Facolta di Ingegneria, Universitd statale di Bergamo, Bergamo, Italy; 
and INFN — Sezione di Milano, Milan, Italy. 

and 

Michel Zamboni-Rached, 

Centro de Ciencias Naturais e Humanas, Universidade Federal do ABC, 

Santo Andre, SP, Brasil 

Abstract — In the First Part of this paper (which is mainly a review) we present sim- 
ple, general and formal, introductions to the ordinary gaussian waves and to the Bessel 
waves, by explicitly separating the case of beams from the case of pulses; and, afterwards, 
an analogous introduction is presented for the Localized Waves (LW), pulses or beams. 
Always we stress the very different characteristics of the gaussian with respect to the 
Bessel waves and to the LWs, showing the numerous important properties of the latter: 
Properties that may find application in all fields in which an essential role is played by 
a wave-equation (like electromagnetism, optics, acoustics, seismology, geophysics, grav- 
itation, elementary particle physics, etc.). The First Part of this review ends with an 
Appendix, wherein: (i) we recall how, in the seventies and eighties, the geometrical meth- 
ods of Special Relativity (SR) predicted — in the sense below specified — the existence of 
the most interesting LWs, i.e., of the X-shaped pulses; and (ii) in connection with the 
circumstance that the X-shaped waves are endowed with Superluminal group-velocities 
(as discussed in the first part of this paper), we briefly mention the various experimen- 
tal sectors of physics in which Superluminal motions seem to appear; in particular, a 
bird's-eye view is presented of the experiments till now performed with evanescent waves 
(and/or tunnelling photons), and with the "localized Superluminal solutions" to the wave 
equations. 

In the Second Part of this work, we address in more detail various theoretical ap- 
proaches leading to nondiffr acting solutions of the linear wave equation in unbounded 
homogeneous media, as well as some interesting applications of these waves. After some 
more introductory remarks (Sec. VI), we analyse in Section VII the general structure of 
the Localized Waves, develop the so called Generalized Bidirectional Decomposition, and 
use it to obtain several luminal and Superluminal nondiffracting solutions of the wave 
equations. In Section VIII we present a method for getting a space-time focusing by a 
continuous superposition of X-Shaped pulses of different velocities. Section IX addresses 
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the properties of chirped optical X-Shaped pulses propagating in material media without 
boundaries. 

Finally, in the Third Part of this paper we "complete" our review by investigating also 
the (not less interesting) case of the suhluminal Localized Solutions to the wave equations, 
which, among the others, will allow us to emphasize the remarkable role of SR, in its 
extended, or rather non-restricted, formulation. [For instance, the various Superluminal 
and subluminal LWs arc expected to be transformed one into the other by suitable Lorentz 
transformations]. We start by studying — by means of various approaches — the very 
peculiar topic of zero-speed waves: Namely, of the localized fields with a static envelope; 
consisting, for instance, in "light at rest" . Actually, in Section X we show how a suitable 
superposition of Bessel beams can be used to construct stationary localized wave fields 
with high transverse localization, and with a longitudinal intensity pattern that assumes 
any desired shape within a chosen interval < 2; < L of the propagation axis. We have 
called Frozen Waves such solutions: As we shall see, they can have a lot of noticeable 
applications. In between, we do not forget to briefly treat the case of not axially- 
symmetric solutions, in terms of higher order Bessel beams. 

In this review we have fixed our attention especially on electromagnetism and optics: 
but results of the present kind are valid, let us repeat, whenever an essential role is 
played by a wave-equation. 

PACS nos.: 03.50.De; 03.30.+P; 03.50.-z; 03.65.Xp; 41.20.Jb; 41.20.-q; 41.85.-p; 
42.25.Bs; 42.25.-p; 42.25.Fx; 43.20.+g; 46.40.-f; 46.40.Cd. 
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FIRST PART 



LOCALIZED WAVES: 
A SCIENTIFIC AND HISTORICAL INTRODUCTION 

1 A General Introduction 

1.1 Preliminary remarks 

Diffraction and dispersion are k;nown since long to be phenomena limiting the applications 
of (optical, for instance) beams or pulses. 

Diffraction is always present, affecting any waves that propagate in two or three- 
dimensional unbounded media, even when homogeneous. Pulses and beams are con- 
stituted by waves travelling along different directions, which produces a gradual spatial 
broadening [Ij. This effect is really a limiting factor whenever a pulse is needed which 
maintains its transverse localization, like, e.g., in free space communications j2j, image 
forming[3\, optical lithography [H [5], electromagnetic iweezers [6], [7j, etcetera. 

Dispersion acts on pulses propagating in material media, causing mainly a temporal 
broadening: An effect known to be due to the variation of the refraction index with the 
frequency, so that each spectral component of the pulse possesses a different phase- velocity. 
This entails a gradual temporal widening, which constitutes a limiting factor when a pulse 
is needed which maintains its time width, like, e.g., in communication systems[H]. 

It is important, therefore, to develop techniques able to reduce those phenomena. 
The so-called localized waves (LW), known also as non-diffracting waves, are indeed able 
to resist diffraction for a long distance in free space. Such solutions to the wave equations 
(and, in particular, to the Maxwell equations, under weak hypotheses) were theoretically 
predicted long time ago[9-12] (cf. also[T3], as well as the Appendix located at the end 
of this First Part), mathematically constructed in more recent times [15] . and soon 
after experimentally produced[16-18]. Today, localized waves are well-established both 
theoretically and experimentally, and are having innovative applications not only in vac- 
uum, but also in material (linear or non-linear) media, showing to be able to resist also 
dispersion. As we were mentioning, their potential applications are being intensively ex- 
plored, always with surprising results, in fields like Acoustics, Microwaves, Optics, and are 
promising also in Mechanics, Geophysics, and even Gravitational Waves and Elementary 
particle physics. Worth noticing appear also the applications of the so-called "Frozen 
Waves" , that will be presented in the Third Part of this work; while rather interesting are 
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the applications already obtained, for instance, in high-resolution ultra-sound scanning of 
moving organs in human body [191 120] . 

To confine ourselves to electromagnetism, let us recall the present-day studies on 
electromagnetic tweezers[21-24], optical (or acoustic) scalpels, optical guiding of atoms 
or (charged or neutral) corpus cles[25-27], optical litography[2Sl EI], optical (or acoustic) 
images^^, communications in free space[30-32,14], remote optical alignment [33]. optical 
acceleration of charged particles, and so on. 

In the following two Subsections we are going to set forth a brief introduction to the 
theory and applications of localized beams and localized pulses, respectively. [34] 

Before going on, let us explicitly remark that — as in any review article, for obvious 
reasons of space — we had to select a few main topics: and such a choice can only be a 
personal one. 

Localized (non- diffracting) beams — The word beam refers to a monochromatic solution 
to the considered wave equation, with a transverse localization of its field. To fix our 
ideas, we shall explicitly refer to the optical case: But our considerations, of course, hold 
for any wave equation (vectorial, spinorial, scalar...: in particular, for the acoustic case 
too). 

The most common type of optical beam is the gaussian one, whose transverse be- 
haviour is described by a gaussian function. But all the common beams suffer a diffraction, 
which spoils the transverse shape of their field, widening it gradually during propagation. 
As an example, the transverse width of a gaussian beam doubles when it travels a distance 
Zdif = v^ttApq/Ao, where Apo is the beam initial width and Aq is its wavelength. One 
can verify that a gaussian beam with an initial transverse aperture of the order of its 
wavelength will already double its width after having travelled a few wavelengths. 

It was generally believed that the only wave devoid of diffraction was the plane wave, 
which does not suffer any transverse changes. Some authors had shown, actually, that it 
isn't the only one. For instance, in 1941 Stratton^lOj obtained a monochromatic solution 
to the wave equation whose transverse shape was concentrated in the vicinity of its prop- 
agation axis and represented by a Bessel function. Such a solution, now called a Bessel 
beam, was not subject to diffraction, since no change in its transverse shape took place 
with time. In ref . |TT] it was later on demonstrated how a large class of equations (including 
the wave equations) admit "non-distorted progressing waves" as solutions; while already 
in 1915, in ref. [12], and subsequently in articles like ref. [35], it was shown the existence of 
soliton-like, wavelet-type solutions to the Maxwell equations. But all such literature did 
not raise the attention it deserved. In the case of ref. [10], this can be partially justified 
since that (Bessel) beam was associated with an infinite power flux [as much as the plane 
waves, incidentally], it being not square-integrable in the transverse direction. An inter- 
esting problem, therefore, was that of investigating what it would happen to the ideal 
Bessel beam solution when truncated by a finite transverse aperture. 

Only in 1987 a heuristical answer came from the known experiment by Durnin et 
al.[36], when it was shown that a realistic Bessel beam, endowed with wavelength Aq = 
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0.6328 /iin and central spolEI Apo = 59 /im, passing through an aperture with radius R = 
3.5 mm is able to travel about 85 cm keeping its transverse intensity shape approximately 
unchanged (in the region p « R surrounding its central peak). In other words, it 
was experimentally shown that the transverse intensity peak, as well as the field in the 
surroundings of it, do not meet any appreciable change in shape all along a large "depth 
of field" . As a comparison, let us recall once more that a gaussian beam with the same 
wavelength, and with the central "spot"[l Apo = 59 /im, when passing through an aperture 
with the same radius R = 3.5 mm doubles its transverse width after 3 cm, and after 6 cm 
its intensity is already diminished by a factor 10. Therefore, in the considered case, a 
Bessel beams can travel, approximately without deformation, a distance 28 times larger 
than a gaussian beam's. 

Such a remarkable property is due to the fact that the transverse intensity fields 
(whose value decreases with increasing p), associated with the rings which constitute the 
(transverse) structure of the Bessel beam, when diffracting end up reconstructing the 
beam itself, all along a large field-depth. This depends on the Bessel beam spectrum 
(wavenumber and frequency) [29| W7\ [33] . as explained in detail in our ref.[38]- Let us 
stress that, given a Bessel and a gaussian beam — both with the same energy E, the same 
spot Apo and passing through apertures with the same radius R in the plane z = — the 
percentage of the total energy E contained inside the central peak region (0 < p < Apo) 
is smaller for a Bessel than for a gaussian beam: This different energy-distribution on 
the transverse plane is responsible for the reconstruction of the Bessel-beam central peak 
even at large distances from the source (and even after an obstacle, provided that its size 
is smaller than the aperture[39-4:l]: a nice property possessed also by the localized pulses 
we are going to examine below [40j). 

It may be worth mentioning that most experiments carried on in this area have been 
performed rapidly and with use, often, of rather simple apparatus: The Durnin et al.'s 
experiment, e.g., had recourse, for the generation of a Bessel beam, to a laser source, an 
annular slit and a lens, as depicted in Fig.([T]). In a sense, such an apparatus produces 
what can be regarded as the cylindrically symmetric generalization of a couple of plane 
waves emitted at angles 6 and —6, with respect to (w.r.t.) the 2;-direction, respectively 
(in which case the plane wave intersection moves along z with the speed c/ cos^). Of 
course, these non-diffracting beams can be generated also by a a conic lens (axicon) [cf., 
e.g., ref. [29j], or by other means like holographic elements [cf., e.g., refs. [33| H2]]. 

Let us stress, as already mentioned at the end of the previous Subsection, that nowa- 
days a lot of interesting applications of non-diffracting beams are being investigated; 
besides the Lu et al.'s ones in Acoustics. In the optical sector, let us recall again those 
of using Bessel beams as optical tweezers able to confine or move around small particles. 
In such theoretical and application noticeable contribution is the one presented 

*Let us define the size of tfie central "spot" of a Bessel beam as the distance, along the transverse 
direction p, at which the first zero occurs of the Bessel function characterizing its transverse shape. 

tin the case of a gaussian beam, let us define the size of its central "spot" as the distance, along the 
transverse direction p, at which its intensity has decayed of the factor 1/e. 
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Figure 1: The simple experimental set-up used by Durnin et al. for generating a Bessel 
beam. 

in refs. [43-45] , wherein, by suitable superpositions of Bessel beams endowed with the 
same frequency but different longitudinal wavenumbers, stationary envelopes have been 
mathematically constructed in closed form, which possess a high transverse localization 
and, more important, a longitudinal intensity-shape that can be freely chosen inside a 
predetermined space-interval < z < L. For instance, a high intensity field, with a 
static envelope, can be created within a tiny region, with negligible intensity elsewhere: 
As already mentioned, the Third Part will deal, among the others, with such "Frozen 
Waves". 

Localized (non-diffracting) pulses — As we have seen in the previous Subsection, the 
existence of non-diffractive (or localized) pulses was predicted since long: cf., once more, 
refs. [121 E] ) and, not less, refs. O |13] , as well as more recent articles like refs. [Ml HT] . 
The modern studies about non-diffractive pulses (to confine ourselves, at least, to the 
ones that attracted more attention) followed a development rather independent of those 
on non-diffracting beams, even if both phenomena are part of the same sector of physics: 
that of Localized Waves. 

In 1983, Brittingham[48j set forth a luminal {V = c) solution to the wave equation 
(more particularly, to the Maxwell equations) which travels rigidly, i.e., without diffrac- 
tion. The solution proposed in ref.[l8] possessed however infinite energy, and once more 
the problem arose of overcoming such a problem. 

A way out was first obtained, as far as we know, by Sezginerpo], who showed how 
to construct finite-energy luminal pulses, which — however — do not propagate without 
distortion for an infinite distance, but, as it is expected, travel with constant speed, and 
approximately without deforming, for a certain (long) depth of field: much longer, in this 
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case too, than that of the ordinary pulses hke the gaussian ones. In a series of subsequent 
papers[30, 31, 50-53], a simple theoretical method was developed, called by those authors 
"bidirectional decomposition", for constructing a new series of non-diffracting luminal 
pulses. 

Eventually, at the beginning of the nineties, Lu et al. [HI [12] constructed, both math- 
ematically and experimentally, new solutions to the wave equation in free space: namely, 
an X-shaped localized pulse, with the form predicted by the so-called extended Special 
Relativity [HI [5l]; for the connection between what Lu et al. called "X-waves" and "ex- 
tended" relativity see, e.g., ref. [I3], while brief excerpts of that theory can be found, for 
instance, in refs. [55-58,15]. Lu et al.'s solutions (which can be called the "classic" ones) 
were continuous superpositions of Bessel beams with the same phase- velocity (i.e., with 
the same axicon angle [59l, [T5l [T^ [5^. alpha); so that they could keep their shape for long 
distances. Such X-shaped waves resulted to be interesting and flexible localized solu- 
tions, and have been afterwards studied in a number of papers, even if their velocity V is 
supersonic or Superluminal (V > c): Actually, when the phase- velocity does not depend 
on the frequency, it is known that such a phase- velocity becomes the "group- velocity"... 
Remembering how a superposition of Bessel beams is generated (for example, by a dis- 
crete or continuous set of annular slits or transducers; or even by a single slit plus a lens), 
it results clear that the energy forming the Localized Waves, coming from those rings, is 
transported at the ordinary speed c of the plane waves in the considered medium[60-62,15] 
(here c, representing the velocity of the plane waves in the medium, is the sound-speed in 
the acoustic case, and the speed of light in the electromagnetic case; and so on). Never- 
theless, the peak of the LWs is faster than c. [Let us explicitly notice that, when using a 
lens after an aperture located at its back focus, as in FiglJl then a classic X-shaped pulse 
can be generated even by a single annular slit, or transducer, illuminated however by a 
light, or sound, pulse: but the previous considerations about the actual transportation- 
speed of the "energy" forming the X-shaped wave remain unaffected. The experimental 
set-up depicted in Fig. 2, with various annular shts, is actually needed only for generating 
(X-shaped, e.g.) pulses more complex than the classic one, namely, depending on the 
co-ordinates z and t not only through the quantity ( = z — Vt: see the following] . 

It is indeed possible to generate (besides the "classic" X-wave produced by Lu et al. in 
1992) infinite sets of new X-shaped waves, with their energy more and more concentrated 
in a spot corresponding to the vertex region |38j. It may therefore appear rather intriguing 
that such a spot [even if no violations of Special Relativity (SR) are obviously implied: all 
the results come from Maxwell equations, or from the wave equations[63l [64]] — travels 
Superluminally when the waves are electromagnetic. For simplicity, we shall call "Super- 
luminal" all the X-shaped waves, even when the waves are acoustic. By Fig. ([3]), which 
refers to an X-wave possessing the velocity > c, we illustrate the fact that, if its vertex 
or central spot is located at Pi at time ti, it will reach the position P2 at a time t + t 
where r = \P2 — Pi\/V < \P2 — Pi\/c: We shall discuss all these points below. 

Soon after having mathematically and experimentally constructed their "classic" 
acoustic X-wave, Lu et al. started applying them to ultrasonic scanning, obtaining — 
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Figure 2: One of the simplest experimental set-ups for generating various kinds of Bessel 
beam superpositions. 




Figure 3: This figure shows an X-shaped wave, that is, a localized Superluminal pulse. 

It refers to an X-wave, possessing the velocity V > c, and illustrates the fact that, if its 
vertex or central spot is located at Pi at time Iq, it will reach the position P2 at a time 
t + T where r = \P2 — Pi\/V < \P2 — Pi\/c: This is something different from the illusory 
"scissor effect" , even if the feeding energy, coming from the regions R, has travelled with 
the ordinary speed c (which is the speed of light in the electromagnetic case, or the sound 
speed in Acoustics, and so on). 
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as we have already said — very high quahty images. Subsequently, in a 1996 e-print 
and report, Recami et al. (see, e.g., ref. [I5] and refs. therein) published the analogous 
X-shaped solutions to the Maxwell equations: By constructing scalar Superluminal local- 
ized solutions for each component of the Hertz potential. That showed, by the way, that 
the localized solutions to the scalar equation can be used, under very weak conditions, 
for obtaining localized solutions to Maxwell's equations too (actually, Ziolkowski et al.|65j 
had found something similar, called by them slingshot pulses, for the simple scalar case; 
but their solution had gone almost unnoticed). In 1997 Saari et al.[17j announced, in 
an important paper, the production in the lab of an X-shaped wave in the optical realm, 
thus proving experimentally the existence of Superluminal electromagnetic pulses. Three 
years later, in 2000, Ranfagni et al. [18] produced, in an experiment of theirs, Superluminal 
X-shaped waves in the microwave region [their paper aroused various criticisms, to which 
those author however responded]. 



2 A More Detailed Introduction 

Let us refer [66] to the differential equation known as homogeneous wave equation: 
simple, but so important in Acoustics, Electromagnetism (Microwaves, Optics,...), Geo- 
physics, and even, as we said, gravitational waves and elementary particle physics: 



92 1 92 \ 

+ + = 0. (1) 



y dx^ dy"^ dz"^ dt"^ ^ 

Let us write it in the cylindrical co-ordinates (p, 0, z) and, for simplicity's sake, confine 
ourselves to axially symmetric solutions ip{p,z;t). Then, eq.([T]) becomes 



92 Id (9^ 1 d'^\ 

+ 77^ + 7^-7^7^ ^(P'^'^) = 0- (2) 



y dp'^ p dp dz"^ (? dt^ 

In free space, solution ip{p,z\ t) can be written in terms of a Bessel- Fourier transform 
w.r.t. the variable p, and two Fourier transforms w.r.t. variables z and t, as follows: 



roo roo roo _ 

ij{p,z,t) = kpJo{kpp)e'''^^e-'^'^{kp,K,uj)dkpdhduj (3) 

JO J —00 J —00 

where Jo(-) is an ordinary zero-order Bessel function and ip{kp, kz,uj) is the transform of 
ilj{p,z,t). 
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Substituting eq.([n]) into eq.([2]), one obtains that the relation, among a;, kp and k^ , 



^ = K + (4) 



has to be satisfied. As a consequence, by using condition (jl]) in eq.Q, any solution to 
the wave equation ([2]) can be written 



^(p,^,t) = r k,Jo{kpp)e'V^^^^^^^'e-'-'S{kp,u)dkpdu; (5) 



where S{kp, uj) is the chosen spectral function, when kz > (and we disregard evenescent 
waves) . 

The general integral solution ([5]) yields for instance the (non-localized) gaussian 
beams and pulses, to which we shall refer for illustrating the differences of the localized 
waves w.r.t. them. 

The Gaussian Beam — A very common (non-localized) beam is the gaussian beam [67]. 
corresponding to the spectrum 



S{kp, uj) = 2a^ e-" 5{u - uq) . (6) 

In eq.(l6]), a is a positive constant, which will be shown to depend on the transverse 
aperture of the initial pulse. 

Figure m illustrates the interpretation of the integral solution ([5]), with spectral func- 
tion ([6]), as a superposition of plane waves. Namely, from FigJUone can easily realize that 
this case corresponds to plane waves propagating in all directions (always with /c^ > 0), the 
most intense ones being those directed along (positive) z. Notice that, in the plane-wave 
case, kz is the longitudinal component of the wave-vector, k = kp + k^, where kp = k^ + ky. 

On substituting eq.([6]) into eq.([5]) and adopting the paraxial approximation, one meets 
the gaussian beam 



2a^ exp ^ 



2(a +iz/2ko) 

where ko = uJq/c. We can verify that such a beam, which suffers transverse diffraction, 
doubles its initial width Apo = 2a after having travelled the distance Zdi{ = \/3koApl/2, 
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Figure 4: Visual interpretation of the integral solution ([5]), with spectral function ([6]), in 
terms of a superposition of plane waves. 

called diffraction length. The more concentrated a gaussian beam happens to be, the 
more rapidly it gets spoiled. 

The Gaussian Pulse — The most common (non-localized) pulse is the gaussian pulse, 
which is got from eq.Q by using the spectrum|68] 




where a and b are positive constants. Indeed, such a pulse is a superposition of gaussian 
beams of different frequency. 

Now, on substituting eq.® into eq.(l5]), and adopting once more the paraxial approx- 
imation, one gets the gaussian pulse: 





(8) 






a 



endowed with speed c and temporal width At = 26, and suffering a progressing enlarge- 
ment of its transverse width, so that its initial value gets doubled already at position 
^^dif = , with Apo = 2a. 



2.1 The localized solutions 



Let us finally go on to the construction of the two most renowned localized waves the 
Bessel beam, and the ordinary X-shaped pulse. 

First of all, it is interesting to observe that, when superposing (axially symmetric) 
solutions of the wave equation in the vacuum, three spectral parameters, (cu, fcp, kz), 
come into the play, which have however to satisfy the constraint (jl]), deriving from the 
wave equation itself. Consequently, only two of them are independent: and we choos^ 
here u and kp. Such a possibility of choosing u) and kp was already apparent in the 
spectral functions generating gaussian beams and pulses, which consisted in the product 
of two functions, one depending only on u and the other on kp. 

We are going to see that further particular relations between u and kp [or, analogously, 
between uj and k^] can be moreover enforced, in order to get interesting and unexpected 
results, such as the localized waves. 

The Bessel beam — Let us start by imposing a linear coupling between uj and kp (it could 
be actually shown [3 7j that it is the unique coupling leading to localized solutions). 
Namely, let us consider the spectral function 



UJ 



S{kp,Lj) 



5{kp sin 6) 

6{u - uJo) 



(10) 



which implies that kp = {u sin 6)/c, with < 6 < 7i/2: A relation that can be regarded as 
a space-time coupling. Let us add that this linear constraint between u and kp, together 
with relation (jl]), yields kz = {ujcos6)/c. This is an important fact, since it has been 
shown elsewhere |66l [38] that an ideal localized wave must contain a coupling of the type 
UJ = Vkz + b, where V and b are arbitrary constants. 

The interpretation of the integral function (I5l), this time with spectrum (ITOl) . as a 
superposition of plane waves is visualized in Figure E) which shows that an axially- 
symmetric Bessel beam is nothing but the result of the superposition of plane waves whose 
wave vectors lay on the surface of a cone having the propagation line as its symmetry axis 
and an opening angle equal to 9; such 9 being called the axicon angle. 

By inserting eq. (fTOl) into eq.(j5]), one gets the mathematical expression of the so-called 
Bessel beam: 



ip{p,z,t) = Jo(^sm9 exp 



I — cos 9 \ z -t 

c \ COS 9 . 



■'■Elsewhere we chose ui and fcr. 
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Figure 5: The axially-symmetric Bessel beam is created by the superposition of plane 
waves whose wave vectors lay on the surface of a cone having the propagation axis as its 
symmetry axis and angle equal to 9 ("axicon angle"). 

This beam possesses phase- velocity t^ph = c/ cos 6', and field transverse shape repre- 
sented by a Bessel function Jo(-) so that its field in concentrated in the surroundings 
of the propagation axis z. Moreover, eq. lfTT]) tells us that the Bessel beam keeps its 
transverse shape (which is therefore invariant) while propagating, with central "spot" 
Ap = 2.405c/ sin ^). 

The ideal Bessel beam, however, is not square-integrable in the transverse direction, 
and is therefore associated with an infinite power fiux: i.e., it cannot be experimentally 
produced. 

But we can have recourse to truncated Bessel beams, generated by finite apertures. 
In this case the (truncated) Bessel beams are still able to travel a long distance while 
maintaining their transfer shape, as well as their speed, approximately unchanged |36 [ |69| 
[70] : That is to say, they still possess a large depth of field. For instance, the field-depth 
of a Bessel beam generated by a circular finite aperture with radius R is given by 



tant^ 

where 9 is the beam axicon angle. In the finite aperture case, the Bessel beam cannot 
be represented any longer by eq. ffTTj) . and one has to calculate it by the scalar diffraction 
theory: By using, for example, Kirchhoff's or Rayleigh-Sommerfeld's diffraction integrals. 
But till the distance .^max one may still use eq. ffTTl) for approximately describing the beam, 
at least in the vicinity of the axis p = 0, that is, for p « R. To realize how much a 
truncated Bessel beam succeeds in resisting diffraction, let us consider also a gaussian 
beam, with the same frequency and central "spot", and compare their field-depths. In 
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particular, let us assume for both beams A = 0.63 jjm and initial central "spot" size 
Apo = 60 yum. The Bessel beam will possess axicon angle 9 = arcsin[2.405c/(ci;Apo)] = 
0.004 rad. Figure [6] depicts the behaviour of the two beams for a circular aperture with 
radius 3.5 mm. We can verify how the gaussian beam doubles its initial transverse 
width already after 3 cm, and after 6 cm its intensity has become an order of magnitude 
smaller. By contrast, the truncated Bessel beam keeps its transverse shape until the 
distance Zmax = R/ tan 6' = 85 cm. Afterwards, the Bessel beam rapidly decays, as a 
consequence of the sharp cut performed on its aperture (such cut being responsible also 
for the intensity oscillations suffered by the beam along its propagation axis, and for the 
fact that eventually the feeding waves, coming from the aperture, at a certain point get 
faint). 




Figure 6: Comparison between a gaussian (a) and a truncated Bessel beam (b). One 
can see that the gaussian beam doubles its initial transverse width already after 3 cm, 
while after 6 cm its intensity decays of a factor 10. By contrast, the Bessel beam does 
approximately keep its transverse shape till the distance 85 cm. 

The zeroth-order (axially symmetric) Bessel beam is nothing but one example of 
localized beam. Further examples are the higher order (not cylindrically symmetric) 
Bessel beams 

ilj(p,(f),z;t) = Ji^i— sin 6 p] exp(iu(f)) exp ( i —cos6 (z ;tM ) 5 (13) 

V c / V c V cosfc* // 

or the Mathieu beams [71], and so on. 

The Ordinary X-shaped Pulse — Following the same procedure adopted in the previous 
subsection, let us construct pulses by using spectral functions of the type 
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5{kp sin 6) 

S{kp,u;) = f F{u) , (14) 

where this time the Dirac delta function furnishes the spectral space-time coupling kp = 
{u!sm9)/c. Function F{u!) is, of course, the frequency spectrum; it is left for the moment 
undetermined. 

On using eq. ffT^ into eq.([n]), one obtains 

ip{p,z,t) = J F{uj) Jq sin 9 exp cos 9 (^z ^^)) ' ^^^^ 

It is easy to see that ip will be a pulse of the type 



^ = tlj{p,z-Vt) (16) 

with a speed V = c/ cos 9 independent of the frequency spectrum F^u). 

Such solutions are known as X-shaped pulses, and are localized (non-diffractive) waves 
in the sense that they do obviously maintain their spatial shape during propagation (see., 
e.g., refs. flS l 138] and refs. therein; as well as the following). 

At this point, some remarkable observations are to be stressed: 

(i) When a pulse consists in a superposition of waves (in this case, Bessel beams) all 
endowed with the same phase- velocity Vph (in this case, with the same axicon angle) 
independent of their frequency, then it is known that the phase-velocity (in this case 
Vph = c/ cos 6') becomes [72] [73] the group- velocity V: That is, V = c/ cos 9 > c. In this 
sense, the X-shaped waves are called "Superluminal localized pulses" (cf., e.g., ref. [T5] 
and refs. therein). [For simplicity, the group- velocity we are talking about [38 | \57 \ [93l [57] 
can be regarded as the peak-velocity. Here, let us only add the observations: (a) that 
the group-velocity for a pulse, in general, is well defined only when the pulse has a clear 
bump in space; but it can be calculated by the approximate, simple relation V ~ du/dk 
only when some extra conditions are satisfied (namely, when u; as a function of k is also 
clearly bumped); and (b) that the group velocity can a priori be evaluated through 
the mentioned, customary derivation of u with respect to the wavenumber for the infinite 
total energy solutions; whilst, for the finite total energy Superluminal solutions, the group- 
velocity cannot be calculated through such an elementary relation, since in those cases it 
does not even exist a one-to-one function uj = uj{kz)]. 
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(ii) Such pulses, even if their group-velocity is Superluminal, do not contradict standard 
physics, having been found in what precedes on the basis of the wave equations — in 
particular, of Maxwell equations [301 US] — only. Indeed, as we shall better see in the 
historical Appendix following below at the end of the First Part, their existence can 
be understood within Special Relativity itself [9, 13, 15, 55-57], on the basis of its ordinary 
Postulates [Mj. Actually, let us repeat it once more, they are fed by waves originating at 
the aperture and carrying energy with the standard speed c of the medium (the light- 
velocity in the electromagnetic case, and the sound-velocity in the acoustic[l6] case). We 
can become convinced about the possibility of realizing Superluminal X-shaped pulses by 
imagining the simple ideal case of a negligibly sized Superluminal source S endowed with 
speed \^ > c in vacuum, and emitting electromagnetic waves W (each one travelling with 
the invariant speed c). The electromagnetic waves will result to be internally tangent to 
an enveloping cone C having S as its vertex, and as its axis the propagation line z of the 
source [511 [T3] : This is completely analogous to what happens for an airplane that moves 
in air with constant supersonic speed. The waves W interfere mainly negatively inside 
the cone C, and constructively on its surface. We can place a plane detector orthogonally 
to z, and record magnitude and direction of the W waves that hit on it, as (cylindrically 
symmetric) functions of position and of time. It will be enough, then, to replace the 
plane detector with a plane antenna which emits — instead of recording — exactly the 
same (axially symmetric) space-time pattern of waves W, for constructing a cone-shaped 
electromagnetic wave C that will propagate with the Superluminal speed V (of course, 
without a source any longer at its vertex...): even if each wave W travels with the invariant 
speed c. Once more, this is exactly what would happen in the case of a supersonic airplane 
(in which case c is the sound speed in air: for simplicity, assume the observer to be at 
rest with respect to the air). For further details, see the quoted references. Actually, by 
suitable superpositions, and interference, of speed-c waves, one can obtain pulses more 
and more localized in the vertex region|38j: That is, very localized field- "blobs" which 
travels with Superluminal group-velocity. This has nothing to do apparently with the 
illusory "scissors effect", since such blobs, along their field-depth, are a priori able, e.g., to 
get two successive (weak) detectors, located at a distance L, to click after a time smaller 
than L/c. Incidentally, an analysis of the above-mentioned case (that of a supersonic 
plane or a Superluminal charge) led, as expected [54j, to the simplest type of "X-shaped 
pulse" [13] . It might be useful, finally, to recall that SR (even the wave-equations have 
an internal relativistic structure!) implies considering also the forward cone: cf. FigjTl 
The truncated X-waves considered in this paper, for instance, must have a leading cone 
in addition to the rear cone; such a leading cone having a role for the peak stability[14]: 
For example, in the approximate case in which we produce a finite conic wave truncated 
both in space and in time, the theory of SR suggested the bi-conic shape (symmetrical in 
space with respect to the vertex 5*) to be a better approximation to a rigidly travelling 
wave (so that SR suggests to have recourse to a dynamic antenna emitting a radiation 
cylindrically symmetric in space and symmetric in time, for a better approximation to an 
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"undistorted progressing wave"). 




Figure 7: The truncated X-waves considered in this paper, as predicted by SR (all wave- 
equations have an intrinsic relativistic structure!), must have a leading cone in addition 
to the rear cone; such a leading cone having a role for the peak stabilityjH]: For example, 
when producing a finite conic wave truncated both in space and in time, the theory of SR 
suggested to have recourse, in the simplest case, to a dynamic antenna emitting a radiation 
cylindrically symmetrical in space and symmetric in time, for a better approximation to 
what Courant and Hilbert[TT] called an "undistorted progressing wave". See the following, 
in the text. 



(iii) Any solutions that depend on z and on t only through the quantity z — Vt, like 
eq. (fT5|) . will appear with a constant shape to an observer travelling along z with the 
speed V. That is, such a solution will propagate rigidly with speed V. This does further 
explain why our X-shaped pulses, after having been produced, will travel almost rigidly 
at speed V (in this case, a faster-than-light group- velocity) , all along their depth of field. 
To be even clearer, let us consider a generic function, depending on z — Vt with V > c, 
and show, by explicit calculations involving the Maxwell equations only, that it obeys the 
scalar wave equation. Following Franco Selleri^Tl], let us consider, e.g., the wave function 



$(x y z t) = (17) 

' ' ' ^J[b - IC{Z - Vt)]^ + (V^ - C2)(x2 + 1/2) 

with a and b non-zero constants, c the ordinary speed of light, and V > c [incidentally, 
this wave function is nothing but the classic X-shaped wave in cartesian co-ordinates]. 
Let us naively verify that it is a solution to the wave equation 



-,9:./ X 1 d'^^(x,y, z,t) , , 

V^$(x, z, t) - - = . (18) 
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On putting 

R = ^[h - ic{z - Vt)f + {V^ - c2)(a;2 + y^) , 
one can write ^ = a/R and evaluate the second derivatives 



(19) 



a d^z i?3 

~ + 3[V ~c ) — 



wherefrom 



and 



1 

a 



i?3 



[b - ic(z - l^t)]^ 

^5 



1 

a 



From the last two equations, remembering the previous definition, one finally gets 



1 



1 



52$ g2^ 

d'^z d^x d'^y (? dH 







that is nothing but the (d'Alembert) wave equation f|T8|) . q.e.d. In conclusion, function 
$ is a solution of the wave equation even if it does obviously represent a pulse (Selleri 
says "a signal"!) propagating with Superluminal speed. 

At this point, the reader should be however warned that all the suhluminal LWs, 
solutions of the ordinary homogeneous wave equation, have appeared till now to present 
sungularities whenever they depend on z and t only via the quantity C, = z — Vt: This 
is still an open, interesting research topic, which is related also to analogous results met 
in gravitation physics. 
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After the previous three important comments, let us go back to our evaluations with 
regard to the X-type solutions to the wave equations. Let us now consider in eq. (fT5l) . for 
instance, the particular frequency spectrum Fioj) given by 



F(a;) = H{u)^ ^^p(-f ^) ' (20) 
where H{ijj) is the Heaviside step-function and a a positive constant. Then, eq. (|T5|) yields 



V^(p,C) = , (21) 

V (« - <f + - 1) 

with C, = z — Vt. This solution fl2T|) is the well-known ordinary, or "classic", X-wave, 
which constitutes a simple example of X-shaped pulse. [TH [15] Notice that function (120|1 
contains mainly low frequencies, so that the classic X-wave is suitable for low frequencies 
only. 

Figure [H] depicts (the real part of) an ordinary X-wave with V = 1.1c and a = 3 m. 




Figure 8: Plot of the real part of the ordinary X-wave, evaluated for V = 1.1c with 
a = 3 m . 

Solutions (fTSjl . and in particular the pulse (12T|) . have got an infinite field-depth, and 
an infinite energy as well. Therefore, as it was done in the Bessel beam case, one should 
pass to truncated pulses, originating from a finite aperture. Afterwards, our truncated 
pulses will keep their spatial shape (and their speed) all along the depth of field 
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z 



R 



(22) 



tan^ ' 



where, as before, R is the aperture radius and 6 the axicon angle. 

Some Further Observations — Let us put forth some further observations. 

It is not strictly correct to call non-diffractive the localized waves, since diffraction af- 
fects, more or less, all waves obeying eq.([l]). However, all localized waves (both beams and 
pulses) possess the remarkable "self-reconstruction" property: That is to say, the localized 
waves, when diffracting during propagation, do immediately re-build their shape[39-41] 
(even after obstacles with size much larger than the characteristic wave-lengths, provided 
it is smaller — as we know — than the aperture size), due to their particular spectral struc- 
ture [as it is shown more in detail, e.g., in the book Localized Waves (J.Wiley; Jan. 2008)]. 
In particular, the "ideal localized waves" (with infinite energy and field-depth) are able 
to re-build themselves for an infinite time; while, as we have seen, the finite-energy (trun- 
cated) ones can do it, and thus resist the diffraction effects, only along a certain depth of 



Let us stress again that the interest of the localized waves (especially from the point 
of view of applications) lies in the circumstance that they are almost non-diffractive, 
rather than in their group-velocity: From this point of view, Superluminal, luminal, and 
subluminal localized solutions are equally interesting and suited to important applications. 

Actually, the localized waves are not restricted to the (X-shaped, Superluminal) ones 
corresponding to the integral solution f|T5|) to the wave equation; and, as we were already 
saying, three classes of localized pulses exist: the Superluminal (with speed V > c), the 
luminal {V = c), and the subluminal {V < c) ones; all of them with, or without, axial 
symmetry, and corresponding in any case to a single, unified mathematical background. 
This issue will be touched again in the present review. 

Incidentally, we have addressed elsewhere topics as: (i) the construction of infinite 
families of generalizations of the classic X-shaped wave [with energy more and more 
concentrated around the vertex: cf., e.g., FigslHl taken from ref.[38]]; as (ii) the behaviour 
of some finite total-energy Superluminal localized solutions (SLS); (iii) the techniques for 
building up new series of SLS's to the Maxwell equations suitable for arbitrary frequencies 
and bandwidths; (iv) questions related with the case of dispersive (and even lossy) 
media; (v) the construction of (infinite or finite energy) Superluminal LWs propagating 
down waveguides or coaxial cables; (vi) finding out Localized Solutions also to equations 
different from the wave equation, as Schroedinger's; (vii) using the above tecniques for 
constructing, in General Relativity, new exact solutions for gravitational ways. In the 
Second Part of this paper we shall come back to some (few) of those points. Let us add 
that X-shaped waves have been easily produced also in nonlinear media [75]. 

A more technical introduction to the subject of localized waves (particularly w.r.t. 
the Superluminal X-shaped ones) can be found in the Second Part of this review, and in 
papers like ref . |55j . 




field... 
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(a) 



(b) 



Figure 9: In Fig. (a) it is represented (in arbitrary units) the square magnitude of the 
"classic", X-shaped Superluminal locahzed solution (SLS) to the wave equation, with 
V = 5c and a = 0.1 m. Families of infinite SLSs however exists, which generalize 
the classic X-shaped solution; for instance, a family of SLSs obtained [38] by suitably 
differentiating the classic X-wave: Fig.(b) depicts the first of them (corresponding to the 
first differentiation) with the same parameters. As we said, the successsive solutions in 
such a family are more and more localized around their vertex. Quantity p is the distance 
in meters from the propagation axis z, while quantity ( is the "V^-cone" variable [ref. [38]] 
(still in meters) ( = z — Vt, with V > c. Since all these solutions depend on z only via the 
variable (, they propagate "rigidly", i.e., as we know, without distortion (and are called 
"localized", or non-diffracting, for such a reason). Here we are assuming propagation in 
the vacuum (or in a homogeneous medium). 

Before going on to the Second Part of this paper, let us end the present First Part by 
a historical (theoretical and experimental) Appendix. 



An APPENDIX to the First Part: 

A HISTORICAL (THEORETICAL AND EXPERIMENTAL) 
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APPENDIX 



In this mainly "historical" Appendix, written as far as possible in a (partially) self- 
consistent form, we shall first refer ourselves, from the theoretical point of view, to the 
most intriguing localized solutions to the wave equation: the Superluminal ones (SLS), 
and in particular the X-shaped pulses. As a start, we shall recall their geometrical in- 
terpretation within SR. Afterwards, to help resolving possible doubts, we shall seize the 
opportunity, given by this Appendix, for presenting a bird's-eye view of the various experi- 
mental sectors of physics in which Superluminal motions seem to appear: In particular, of 
the experiments with evanescent waves (and/or tunnelling photons), and with the SLS's 
we are more interested in here. In some parts of this Appendix the propagation-line is 
called X, and no longer z, without originating, however, any interpretation problems. 



3 An Introduction to the APPENDIX 

The question of Superluminal {V"^ > c^) objects or waves has a long story. Still in pre- 
relativistic times, one meets various relevant papers, from those by J.J.Thomson to the 
interesting ones by A.Sommerfeld. It is well-known, however, that with SR the conviction 
spread out that the speed c of light in vacuum was the upper limit of any possible speed. 
For instance, R.C.Tolman in 1917 believed to have shown by his "paradox" that the 
existence of particles endowed with speeds larger than c would have allowed sending 
information into the past. Our problem started to be tackled again only in the fifties 
and sixties, in particular after the papers[76] by E.C.George Sudarshan et al., and, later 
on [771 [78] . by one of the present authors with R.Mignani et al., as well as — to confine 
ourselves at present to the theoretical researches — by H.C.Corben and others. The first 
experimental attempts were performed by T.Alvager et al. 

We wish to face the still unusual issue of the possible existence of Superluminal 
wavelets, and objects — within standard physics and SR, as we said — since at least four 
different experimental sectors of physics seem to support such a possibility [apparently 
confirming some long-standing theoretical predictions [541 [9l [761 EH]]- The experimental 
review will be necessarily short, but we shall provide the reader with enough bibliograph- 
ical information, limited for brevity's sake to the last century only (i.e., up-dated till the 
year 2000 only). 
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4 APPENDIX: Historical Recollections - Theory 

A simple theoretical framework was long ago proposed [76| [5^ [77], merely based on 
the space-time geometrical methods of SR, which appears to incorporate Superluminal 
waves and objects, and in a sense predicts[9] among the others the Superluminal X-shaped 
waves, without violating the Relativity principles. A suitable choice of the Postulates of 
SR (equivalent of course to the other, more common, choices) is the following one: (i) the 
standard Principle of Relativity; and (ii) space-time homogeneity and space isotropy. It 
follows that one and only one invariant speed exists; and experience shows that invariant 
speed to be the light-speed, c, in vacuum: The essential role of c in SR being just due 
to its invariance, and not to the fact that it be a maximal, or minimal, speed. No sub- 
or Super-luminal objects or pulses can be endowed with an invariant speed: so that their 
speed cannot play in SR the same essential role played the light-speed c in vacuum. 
Indeed, the speed c turns out to be also a limiting speed: but any limit possesses two 
sides, and can be approached a priori both from below and from above: See FiglTOl As 
E.C.G.Sudarshan put it, from the fact that no one could climb over the Himalayas ranges, 
people of India cannot conclude that there are no people North of the Himalayas... Indeed, 
speed-c photons exist, which are born, live and die just "at the top of the mountain," 
without any need for performing the impossible task of accelerating from rest to the 
light-speed. [Actually, the ordinary formulation of SR has been too much restricted: For 
instance, even leaving Superluminal speeds aside, it can be easily so widened as to include 
antimatter [511 EH ET]]. 




Figure 10: Energy of a free object as a function of its speed. [76l EZl |M] 

An immediate consequence is that the quadratic form c^dt^ — da;^ = dx^dx'^, called 
ds^, with da?^ = dx^ + dy^ + dz^, results to be invariant, except for its sign. Quantity 
ds^, let us recall, is the four-dimensional length-element square, along the space-time 
path of any object. In correspondence with the positive (negative) sign, one gets the 
subluminal (Superluminal) Lorentz "transformations" [LT] . More specifically, the ordinary 
subluminal LTs are known to leave, e.g., the quadratic forms dx^dx^, dp^dp^ and dx^dp'^ 
exactly invariant, where the are the component of the energy-impulse four- vector; 
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while the Superluminal LTs, by contrast, have to change (only) the sign of such quadratic 
forms. This is enough for deducing some important consequences, like the one that a 
Superluminal charge has to behave as a magnetic monopole, in the sense specified in 
ref.[5lj and refs. therein. 

A more important consequence, for us, is — see FigHT] — that the simplest subluminal 
object, namely a spherical particle at rest (which appears as ellipsoidal, due to Lorentz 
contraction, at subluminal speeds v), will appear P [SU [15] as occupying the cylindrically 
symmetrical region bounded by a two-sheeted rotation hyperboloid and an indefinite 
double cone, as in Fig. (lid), for Superluminal speeds V. In the limiting case of a point- 
like particle, one obtains only a double cone. 




(a) (b) (c) (d) 

Figure 11: An intrinsically spherical (or pointlike, at the limit) object appears in the 
vacuum as an ellipsoid contracted along the motion direction when endowed with a speed 
V < c. By contrast, if endowed with a speed V > c (even if the c-speed barrier cannot 
be crossed, neither from the left nor from the right), it would appear P, El] no longer as 
a particle, but rather as an "X-shaped" wave travelling rigidly: Namely, as occupying 
the region delimited by a double cone and a two-sheeted hyperboloid — or as a double 
cone, at the limit-, and moving without distortion in the vacuum, or in a homogeneous 
medium, with Superluminal speed V [the square cotangent of the cone semi-angle being 
{V/cy — 1]. For simplicity, a space axis is skipped. This figure is taken from refs. [H El] . 

Such a result is got by writing down the equation of the world-tube of a subluminal parti- 
cle, and transforming it simply by changing the sign of the quadratic forms entering that 
equation. Thus, in 1980-1982, it was predicted[9] that the simplest Superluminal object 
appears (not as a particle, but as a field or rather) as a wave: namely, as an "X-shaped 
pulse", the cone semi-angle a being given (with c = 1) by cotga = ^/V^'^^. Such 
X-shaped pulses will move rigidly with speed V along their motion direction: In fact, any 
"X-pulse" can be regarded at each instant of time as the (Superluminal) Lorentz transform 
of a spherical object, which of course moves in vacuum — or in a homogeneous medium — 
without any deformation as time elapses. The three-dimensional picture of Fig. (lid) 
appears in FigIT2l where its annular intersections with a transverse plane are shown (cf. 
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refs.[n]). The X-shaped waves here considered are merely the simplest ones: if one starts 
not from an intrinsically spherical or point-like object, but from a non-spherically sym- 
metric particle, or from a pulsating (contracting and dilating) sphere, or from a particle 
oscillating back and forth along the motion direction, then their Superluminal Lorentz 
transforms would result to be more and more complicated. The above-seen "X-waves", 
however, are typical for a Superluminal object, as much as the spherical or point-like 
shape is typical, let us repeat, for a subluminal object. 




Figure 12: Here we show the intersections of the Superluminal object T represented in 
Fig. (lid) with planes P orthogonal to its motion line (the z-axis). For simplicity, we 
assumed again the object to be spherical in its rest-frame, and the cone vertex C to 
coincide with the origin O for t = 0. Such intersections evolve in time so that the same 
pattern appears on a second plane — shifted by Ax — after the time At = Ax/V. On 
each plane, as time elapses, the intersection is therefore predicted by (extended) SR to 
be a circular ring which, for negative times, goes on shrinking until it reduces to a circle 
and then to a point (for t = 0); afterwards, such a point becomes again a circle and then 
a circular ring that goes on broadening[9l EH [15]. This picture is taken from refs. [HI El] • 
[Notice that, if the object is not spherical when at rest (but, e.g., is ellipsoidal in its own 
rest-frame), then the axis of T will no longer coincide with x, but its direction will depend 
on the speed V of the tachyon itself]. For the case in which the space extension of the 
Superluminal object T is finite, see refs.jH] 

Incidentally, it has been believed for a long time that Superluminal objects would 
have allowed sending information into the past; but such problems with causality seem to 
be solvable within SR. Once SR is generalized in order to include Superluminal objects 
or pulses, no signal travelling backward in time is apparently left. For a solution of those 
causal paradoxes, see refs. [581 EH [ZS] and references therein. 

When addressing the problem, within this elementary Appendix, of the production of 
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an X-shaped pulse like the one depicted in Fig{T2] (maybe truncated, in space and in time, 
by use of a finite antenna radiating for a finite time), all the considerations expounded 
under point (ii) of the subsection The Ordinary X-shaped Pulse become in order: And, 
here, we simply refer to them. Those considerations, together with the present ones 
(related, e.g., to FigJT^. suggest the simplest antenna to consist in a series of concentric 
annular slits, or transducers [like in Figj2], which suitably radiate following specific time 
patterns: See, e.g., refs.[79] and refs. therein. Incidentally, the above procedure can lead 
to a very simple type of X-shaped wave, as investigated below. 
From the present point of view, it is rather interesting to note that, during the last 
fifteen years, X-shaped waves have been actually found as solutions to the Maxwell and 
to the wave equations [let us repeat that the form of any wave equations is intrinsically 
relativistic] . In order to see more deeply the connection existing between what predicted by 
SR (see, e.g., Figs JTT1[T2!) and the localized X-waves mathematically, and experimentally, 
constructed in recent times, let us tackle below, in detail, the problem of the (X-shaped) 
field created by a Superluminal electric charge [T3]. by following a paper recently appeared 
in Physical Review eH 



4.1 The particular X-shaped field associated with a Superlumi- 
nal charge 

It is well-known by now that Maxwell equations admit of wavelet-type solutions endowed 
with arbitrary group- velocities (0 < < oo). We shall again confine ourselves, as above, 
to the localized solutions, rigidly moving: and, more in particular, to the Superluminal 
ones (SLS), the most interesting of which resulted to be X-shaped, as we have already 
seen. The SLSs have been actually produced in a number of experiments, always by 
suitable interference of ordinary-speed waves. In this subsection we show, by contrast, 
that even a Superluminal charge creates an electromagnetic X-shaped wave, in agreement 
with what predicted [9l [Mj within SR. In fact, on the basis of Maxwell equations, one is 
able to evaluate the field associated with a Superluminal charge (at least, under the rough 
approximation of pointlikeness) : As announced in what precedes, it results to constitute 
a very simple example of true X-wave. 

Indeed, the theory of SR, when based on the ordinary Postulates but not restricted 
to subluminal waves and objects, i.e., in its extended version, predicted the simplest 
X-shaped wave to be the one corresponding to the electromagnetic field created by a 
Superluminal chargejSOl [13]. It seems really important evaluating such a field, at least 
approximately, by following ref . [I3] . 

§At variance with the old times — e.g., at the beginning of the seventies our papers on similar subjects 
were always rejected by the most important journals — , things have now changed as to superluminal 
motions: For instance, the paper of ours quoted in Ref. [13], submitted in 2002 to PRL, was diverted to 
PRE, but was eventually published therein in 2004, even if dealing — as we said — with the X-shaped field 
generated by a superhuninal electric charge... 



26 



The toy-model of a pointlike Superluminal charge — Let us start by considering, formally, 
a pointlike Superluminal charge, even if the hypothesis of pointlikeness (already unaccept- 
able in the sub luminal case) is totally inadequate in the Superluminal case|54j- Then, 
let us consider the ordinary vector-potential A^^ and a current density = (0,0, j^;^") 
flowing in the 2;-direction (notice that the motion line is still the axis z). On assuming 
the fields to be generated by the sources only, one has that = (O,O,A^;0), which, 
when adopting the Lorentz gauge, obeys the equation = j^. We can write such non- 
homogeneous wave equation in the cylindrical co-ordinates (p, 9, z] t); for axial symmetry 
[which requires a priori that A^ = A^{p, z]t)], when choosing the 'V-cone variables" 
( = z — Vt; T] = z + Vt , with > c^, we arrive [T3] at the equation 



-P 



d ( d 
P 



+ 



1 52 



+ 



1 92 



dp \ dp I 7^ dC^ 7^ drj'^ d(dri 



A^{pX,v) = fipX,v) 



(23) 



where p assumes the two values p = 3,0 only, so that = (0,0,^42;^), and 7^ = 
[V'^ — [Notice that, whenever convenient, we set c = 1]. Let us now suppose A^ 

to be actually independent of rj, namely, A'^ = A'^(p, C)- Due to eq. (!23|) . we shall have 



= i^{p, C) too; and therefore jz = Vj° (from the continuity equation), and A^ = V(j)/c 
(from the Lorentz gauge). Then, by calling if) = A^, we end up in two equations |13j. 
which allow us to analyse the possibility and consequences of having a Superluminal 
pointlike charge, e, travelling with constant speed V along the z-axis (p = 0) in the 
positive direction, in which case jz = 6{p) / p 5(C)- Indeed, one of those two equations 
becomes the hyperbolic equation 
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which can be solved [I3] in few steps. First, by applying (with respect to the variable p) 

POD 

the Fourier-Bessel (FB) transformation fix) = / f2/(f2) Jo(fix) dfi, quantity Jo(fix) 

JO 

being the ordinary zero-order Bessel function. Second, by applying the ordinary Fourier 
transformation with respect to the variable ( (going on, from (, to the variable u). And, 
third, by finally performing the corresponding inverse Fourier and FB transformations. 
Afterwards, it is enough to have recourse to formulae (3.723.9) and (6.671.7) of ref. [81]. 
still with ( = z — Vt , for being able to write down the solution of eq. (1241) in the form 
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In FigHSlwe show our solution = ip, as a. function of p and (, evaluated for 7 = 1 (i.e., 
for V = cV2). Of course, we skipped the points at which must diverge, namely the 
vertex and the cone surface. 




-0.4 



-0.6 



-0.8 



p(nm) 



-1 



Figure 13: Behaviour of the field = generated by a charge supposed to be Super- 
luminal, as a function of p and = z — Vt, evaluated for 7 = 1 (i.e., for V = cV2y. 
According to ref.|T3] [Of course, we skipped the points at which ip must diverge: namely, 
the vertex and the cone surface]. 

For comparison, one may recall that the classic X-shaped solution [H] of the homoge- 
neous wave-equation — which is shown, e.g., in FigsIHl El [12] — has the form (with a > 0): 



The second one of eqs.(25) includes expression (26), given by the spectral parameter [38l [82] 
a = 0, which indeed corresponds to the non-homogeneous case [the not negligible fact that 
for a = these equations differ for an imaginary unit [SU [83] will be discussed elsewhere]. 

It is rather important, at this point, to notice that such a solution, eq.(25), does 
represent a wave existing only inside the (unlimited) double cone C generated by the 
rotation around the 2;-axis of the straight lines p = ±7^: This too is in full agreement 
with the predictions of the extended theory of SR. For the explicit evaluation of the 
electromagnetic fields generated by the Superluminal charge (and of their boundary values 
and conditions) we confine ourselves here to merely quoting ref.[l3]. Incidentally, the 
same results found by following the above procedure can be obtained by starting from the 
four-potential associated with a subluminal charge (e.g., an electric charge at rest), and 
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then applying to it the suitable Superluminal Lorentz "transformation" . One should also 
notice that this double cone does not have much to do with the Cherenkov cone [HH EDI El]; 
and that a Superluminal charge travelling at constant speed, in the vacuum, does not lose 
energy: See, e.g., FigHH [which reproduces figure 27 at page 80 of ref.jM]]- 
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Figure 14: The spherical equipotential surfaces of the electrostatic field created by a 
charge at rest get transformed into two-sheeted rotation-hyperboloids, contained inside 
an unlimited double-cone, when the charge travels at Superluminal speed (cf. refs.plSlIM]). 
This figures shows, among the others, that a Superluminal charge travelling at constant 
speed, in a homogeneous medium like the vacuum, does not lose energy [80j. Let us 
mention, incidentally, that this double cone has nothing to do with the Cherenkov cone [84] . 
[The present picture is a reproduction of figure 27, appeared in 1986 at page 80 of ref.[M]]. 

Outside the cone C, i.e., for < 7 | C |< we get as expected no field, so that 
one meets a field discontinuity when crossing the double-cone surface. Nevertheless, the 
boundary conditions imposed by Maxwell equations are satisfied by our solution (25), since 
at each point of the cone surface the electric and the magnetic field are both tangent to 
the cone: also for a discussion of this point we refer to quotation [13j. 

Here, let us stress that, when \^ — > cxd, and therefore 7 — >• 0, the electric field tends 
to vanish, while the magnetic field tends to the value = — vre/p^: This does agree 
once more with what expected from extended SR, which predicted Superluminal charges 
to behave (in a sense) as magnetic monopoles. In the present paper we can only mention 
such a circumstance, and refer to citations[77l [511 EHl [83], and papers quoted therein. 
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5 APPENDIX: A Glance at the Experimental State- 
of-the-Art 



Extended relativity can allow a better understanding of many aspects also of ordi- 
nary physics[5l], even if Superluminal objects (tachyons) did not exist in our cosmos as 
asymptotically free objects. Anyway, at least three or four different experimental sectors 
of physics seem to suggest the possible existence of faster-than-light motions, or, at least, 
of Superluminal group-velocities. We are going to put forth in the following some infor- 
mation about the experimental results obtained in two of those different physics sectors, 
with a mere mention of the others. 

Neutrinos - First: A long series of experiments, started in 1971, seems to show that the 
square mo^ of the mass mo of muon-neutrinos, and more recently of electron-neutrinos too, 
is negative; which, if confirmed, would mean that (when using a naive language, commonly 
adopted) such neutrinos possess an "imaginary mass" and are therefore tachyonic, or 
mainly tachyonic|85 [ [5H [86] . [In extended SR, however, the dispersion relation for a free 
Superluminal object does become u"^ — k'^ = — or — = —m"^, and there is no 
need at all, therefore, of imaginary masses]. 

Galactic Micro-quasars - Second: As to the apparent Superluminal expansions ob- 
served in the core of quasars[87] and, recently, in the so-called galactic micro- quasars [88], 
we shall not really deal with that problem, too far from the other topics of this paper; 
without mentioning that for those astronomical observations there exist also orthodox 
interpretations, based on ref. [89], that are still accepted by the majority of the astro- 
physicists. For a theoretical discussion, see ref. [90]. Here, let us only emphasize that 
simple geometrical considerations in Minkowski space show that a single Superluminal 
source of light would appear [901 El]- (i) initially, in the "optical boom" phase (analo- 
gous to the acoustic "boom" produced by an airplane traveUing with constant supersonic 
speed), as an intense source which suddenly comes into view; and which, afterwards, (ii) 
seems to split into TWO objects receding one from the other with speed V > 2c [all this 
being similar to what has been actually observed, according to refs.|88j]. 

Evanescent waves and "tunnelling photons" - Third: Within quantum mechanics (and 
precisely in the tunnelling processes), it had been shown that the tunnelling time — firstly 
evaluated as a simple Wigner's "phase time" and later on calculated through the analysis 
of the wavepacket behaviour — does not depend [9T| 192] on the barrier width in the case 
of opaque barriers ("Hartman effect"). This implies Superluminal and arbitrarily large 
group-velocities V inside long enough barriers: see FiglTSl 

Experiments that may verify this prediction by, say, electrons or neutrons are difficult 
and rare [93l m] . Luckily enough, however, the Schroedinger equation in the presence of 
a potential barrier is mathematically identical to the Helmholtz equation for an electro- 
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Figure 15: Behaviour of the average "penetration time" (in seconds) spent by a tunnelhng 
wavepacket, as a function of the penetration depth (in angstroms) down a potential barrier 
(from Olkhovsky et al., ref.[92]). According to the predictions of quantum mechanics, 
the wavepacket speed inside the barrier increases in an unhmited way for opaque barriers; 
and the total tunnelling time does not depend on the barrier width. pTl W2\ 

magnetic wave propagating, for instance, down a metallic waveguide (along the 2;-axis): 
as shown, e.g., in refs.[95]; and a barrier height U bigger than the electron energy E 
corresponds (for a given wave frequency) to a waveguide of transverse size lower than 
a cut-off value. A segment of "undersized" guide — to go on with our example — does 
therefore behave as a barrier for the wave {photonic barrier), as well as any other photonic 
band-gap filters. The wave assumes therein — like a particle inside a quantum barrier — 
an imaginary momentum or wavenumber and, as a consequence, results exponentially 
damped along x [see, e.g. FigHE]: It becomes an evanescent wave (going back to normal 
propagation, even if with reduced amplitude, when the narrowing ends and the guide 
returns to its initial transverse size). Thus, a tunnelling experiment can be simulated 
by having recourse to evanescent waves (for which the concept of group velocity can be 
properly extended: see the first one of refs.[57j). 

The fact that evanescent waves travel with Superluminal speeds (cf., e.g., Fig JTTI) has 
been actually verified in a series of famous experiments. Namely, various experiments, 
performed since 1992 onwards by G.Nimtz et al. in Cologne [96]. by R.Chiao, P.G.Kwiat 
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Figure 16: Picture of the damping taking place inside a barrier (from ref.[57]): this 
damping does reduce the amphtude of the tunnelhng wavepacket, imposing a practical 
limit on the adoptable barrier length. 

and A.Steinberg at Berkeley [97], by A.Ranfagni and colleagues in Florence[T8], and by 
others in Vienna, Orsay, Rennes, etcetera[98j, verified that "tunnelling photons" travel 
with Superluminal group velocities [Such experiments raised a great deal of interest [99]. 
also within the non-specialized press, and were reported in Scientific American, Nature, 
New Scientist, etc.]. Let us further remark that also extended SR had predicted |100j 
evanescent waves to be endowed with faster-than-c speeds; the whole matter appears to 
be therefore theoretically self consistent. The debate in the current literature does not 
refer to the experimental results (which can be correctly reproduced even by numerical 
simulations [63l [64] based on Maxwell equations only: Cf. Figs IT8lfT9l) . but rather to 
the question whether they allow, or do not allow, sending signals or information with 
Superluminal speed (see, e.g., refs. [101j ). 

In the above-mentioned experiments one meets a substantial attenuation of the con- 
sidered pulses — cf. FigJlH] — during tunnelling (or during propagation in an absorbing 
medium) : However, by employing "gain doublets" , it has been recently reported the ob- 
servation of undistorted pulses propagating with Superluminal group- velocity with a small 
change in amplitude (see, e.g., ref. |103] ). 

Let us emphasize that some of the most interesting experiments of this series seem 
to be the ones with TWO or more "barriers" (e.g., with two gratings in an optical 
fiber [102] . or with two segments of undersized waveguide separated by a piece of normal- 
sized waveguide ^104j : Fig j20|) . 

For suitable frequency bands — namely, for "tunnelling" far from resonances — , it was 
found by us that the total crossing time does not depend on the length of the intermedi- 
ate (normal) guide: that is, that the beam speed along it is infinite [105[ 11041 [93] . This 
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Figure 17: Simulation of tunnelling by experiments with evanescent classical waves (see 
the text), which were predicted to be Superluminal also on the basis of extended SR |100] . 
The figure shows one of the measurement results by Nimtz et al.[96]; that is, the average 
beam speed while crossing the evanescent region ( = segment of undersized waveguide, or 
"barrier") as a function of its length. As theoretically predicted [OTl 1100] . such an average 
speed exceeds c for long enough "barriers". Further results appeared in ref. [102j . and are 
reported below: see FigsJSOl and [2T] in the following . 

does agree with what predicted by Quantum Mechanics for the non-resonant tunnelling 
through two successive opaque barriers [105j : Figj2Tl Such a prediction has been verified 
first theoretically, by Y.Aharonov et al. |105j . and then, a second time, experimentally: by 
taking advantage of the circumstance that evanescence regions can consist in a variety of 
photonic band-gap materials or gratings (from multilayer dielectric mirrors, or semicon- 
ductors, to photonic crystals). Indeed, the best experimental confirmation has come by 
having recourse to two gratings in an optical fiber |102] : see Figsj22]and 23; in particular, 
the rather peculiar (and quite interesting) results represented by the latter. 
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Figure 18: The delay of a wavepacket crossing a barrier (cf., e.g., FiglT7]is due to the initial 
discontinuity. We then performed suitable numerical simulations [63] by considering an 
(indefinite) undersized waveguide, and therefore eliminating any geometric discontinuity 
in its cross-section. This figure shows the envelope of the initial signal. Inset (a) depicts 
in detail the initial part of this signal as a function of time, while inset (b) depicts the 
gaussian pulse peak centered at t = 100 ns . 
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Figure 19: Envelope of the signal in the previous figure (Fig JTSl) after having travelled a 
distance L = 32.96 mm through the mentioned undersized waveguide. Inset (a) shows 
in detail the initial part (in time) of such arriving signal, while inset (b) shows the peak 
of the gaussian pulse that had been initially modulated by centering it at t = 100 ns. 
One can see that its propagation took zero time, so that the signal travelled with infinite 
speed. The numerical simulation has been based on Maxwell equations only. Going on 
from Fig. 18 to Fig. 19 one verifies that the signal strongly lowered its amplitute: However, 
the width of each peak did not change (and this might have some relevance when thinking 
of a Morse alphabet "transmission": see the text) . 



We cannot skip a further topic — which, being delicate, should not appear, probably, in 
a brief overview like this — since it is presently arising more and more interest [103] . Even 
if all the ordinary causal paradoxes seem to be solvable [58l [5^ |57] . nevertheless one has to 
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Figure 20: Very interesting experiments have been performed with TWO successive barri- 
ers, i.e., with two evanescence regions: For example, with two gratings in an optical fiber. 
This figure [57j refers to the interesting experiment |104] performed with microwaves trav- 
elling along a metallic waveguide: the waveguide being endowed with two classical barriers 
(undersized guide segments). See the text . 
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Figure 21: Scheme of the (non- resonant) tunnelling process, through two successive 
(opaque) quantum barriers. Far from resonances, the (total) phase time for tunnelling 
through the two potential barriers does depend neither on the barrier widths nor on the 
distance between the barriers ("generalized Hartman effect" ) [T05| [93 | [T06] . See the text. 
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Figure 22: Realization of the quantum-theoretical set-up represented in FigllH by using, 
as classical (photonic) barriers, two gratings in an optical fiber |106] . The corresponding 
experiment has been performed by Longhi et al. |102] 



bear in mind that (whenever it is met an object, (9, travelling with Superluminal speed) 
one may have to deal with negative contributions to the tunnelling times p^Tf [SU [93] : 
and this should not be regarded as unphysical. In fact, whenever an "object" (particle, 
electromagnetic pulse,,...) O ot'ercomes [581 EH the infinite speed with respect to a certain 
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Figure 23: Off-resonance tunnelling time versus barrier separation for the rectangular 
symmetric DB FBG structure considered in ref.[102j (see Fig.22). The solid line is the 
theoretical prediction based on group delay calculations; the dots are the experimental 
points as obtained by time delay measurements [the dashed curve is the expected tran- 
sit time from input to output planes for a pulse tuned far away from the stopband of 
the FBGs]. The experimental results p^2] do confirm — as well as the early ones in 
refs.[TD4j] — the theoretical prediction of a "generalized Hartman Effect": in particular, 
the independence of the total tunnelling time from the distance between the two barriers. 



observer, it will afterwards appear to the same observer as the "aniz-object" O travelling 
in the opposite space direction [76| [Ml [58] . For instance, when going on from the lab 
to a frame moving in the same direction as the particles or waves entering the barrier 
region, the object O penetrating through the final part of the barrier (with almost infinite 
speed [92], [911 [Ml [93], like in Figs. 15) will appear in the frame JF as an anti-object O crossing 
that portion of the barrier in the opposite space- direction^^ [Ml [ZS]- In the new frame 
JF, therefore, such anti-object O would yield a negative contribution to the tunnelling 
time: which could even result, in total, to be negative. For any clarifications, see the 
quoted references. Let us stress, here, that even the appearance of such negative times 
had been predicted within SR itself |107] . on the basis of its ordinary postulates; and has 
been recently confirmed by quantum-theoretical evaluations too [93l [T08] . (In the case of 
a non-polarized beam,, the wave anti-packet coincides with the initial wave packet; if a 
photon is however endowed with helicity A = +1, the anti-photon will bear the opposite 
helicity A = —1). From the theoretical point of view, besides the above-quoted papers 
(in particular refs. [93l [9T]). see more specifically refs. |109j . On the (very interesting!) 
experimental side, see the intriguing papers [llOJ. 

Let us add here that, via quantum interference effects, it is possible to obtain di- 
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electrics with refraction indices very rapidly varying as a function of frequency, also in 
three-level atomic systems, with almost complete absence of light absorption (i.e., with 
quantum induced transparency) The group velocity of a light pulse propagating in 
such a medium can decrease to very low values, either positive or negative, with no pulse 
distortion. It is known that experiments have been performed both in atomic samples 
at room temperature, and in Bose-Einstein condensates, which showed the possibility of 
reducing the speed of light to a few meters per second. Similar, but negative group 
velocities, implying a propagation with Superluminal speeds thousands of time higher 
than the previously mentioned ones, have been recently predicted also in the presence 
of such an "electromagnetically induced transparency", for light moving in a rubidium 
condensate [112]. Finally, let us recall that faster-than-c propagation of light pulses can 
be (and has been, in same cases) observed also by taking advantage of the anomalous dis- 
persion near an absorbing line, or nonlinear and linear gain lines — as already seen — , or 
nondispersive dielectric media, or inverted two-level media, as well as of some parametric 
processes in nonlinear optics (cf., e.g., G.Kurizki et al.'s works). 

D) Superluminal Localized Solutions (SLS) to the wave equations. The "X- 
shaped waves" - The fourth sector (to leave aside the others) is not less important. It 
came into fashion again, when it was rediscovered in a series of remarkable works that 
any wave equation — to fix the ideas, let us think of the electromagnetic case — admit also 
solutions as much sub-luminal as Super-luminal (besides the luminal ones, having speed 
c/n). Let us recall, indeed, that, starting from pioneering works as H.Bateman's, it had 
slowly become known that all wave equations admit soliton-like (or rather wavelet-type) 
solutions with sub-luminal group velocities. Subsequently, also Superluminal solutions 
started to be written down (in one case [35] just by the mere application of a Superluminal 
Lorentz "transformation" [54]). 

As we know, a remarkable feature of some new solutions of these (which attracted 
much attention for their possible applications) is that they propagate as localized, non- 
diffracting pulses, also because of their self-reconstruction property. It is easy to realize 
the practical importance, for instance, of a radio transmission carried out by localized 
beams, independently of their speed; but non-diffracting wave packets can be of use even 
in theoretical physics for a reasonable representation of elementary particles; and so on. 
Incidentally, from the point of view of elementary particles, it can be a source of meditation 
the fact that the wave equations possess pulse- type solutions that, in the sub luminal case, 
are ball-like (cf. Fig l24l) : this can have a bearing on the corpuscle/wave duality problem 
met in quantum physics (besides agreeing, e.g., with Fig JTTl) . Further comments on this 
point are to be found below. 

At the cost of repeating ourselves, let us emphasize once more that, within extended 
SR, since 1980 it had been found that — whilst the simplest subluminal object conceivable 
is a small sphere, or a point in the limiting case — the simplest Superluminal objects 
results by contrast to be (see refs.[9], and FigsJTTl and fT2] of this paper) an "X-shaped" 
wave, or a double cone as its limit, which moreover travels without deforming — i.e.. 
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Figure 24: The wave equations possess pulse-type solutions that, in the subluminal case, 
are ball-like, in agreement with FigJTTl For comments, see the text. 

rigidly — in a homogeneous medium. It is not without meaning that the most interesting 
localized solutions to the wave equations happened to be just the Superluminal ones, 
and with a shape of that kind. Even more, since from Maxwell equations under simple 
hypotheses one goes on to the usual scalar wave equation for each electric or magnetic field 
component, one expected the same solutions to exist also in the field of acoustic waves, 
of seismic waves, and of gravitational waves too: and this has already been demonstrated 
in the literature for all those cases, and especially in Acoustics. Actually, such pulses 
(as suitable superpositions of Bessel beams) were mathematically constructed for the first 
time, by Lu et al. in Acoustics: and were then called "X-waves" or rather X-shaped 
waves. (One should not forget that, however, LWs have been constructed in exact form 
even for other equations, as Schroedinger's and as Einstein's). 

It is indeed important for us that the X-shaped waves have been really produced in 
experiments, both with acoustic and with electromagnetic waves; that is, X-pulses were 
produced that, in their medium, travel undistorted with a speed larger than sound, in 
the first case, and than light, in the second case. In Acoustics, the first experiment was 
performed by Lu et al. themselves in 1992, at the Mayo Clinic (and their papers received 
the first IEEE 1992 award). In the electromagnetic case, certainly more intriguing, 
Superluminal localized X-shaped solutions were first mathematically constructed (cf., e.g., 
Fig j25l) in refs.|15j. and later on experimentally produced by Saari et al.[T7] in 1997 
at Tartu by visible light (Fig j26|l . and more recently by Ranfagni et al. at Florence 
by microwaves fT8]. In the theoretical sector the activity has been not less intense, in 
order to build up — for example — analogous new solutions with finite total energy or 
more suitable for high frequencies, on one hand, and localized solutions Superluminally 
propagating even along a normal waveguide (cf. Fig|5]), on another hand, and so on. 
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Figure 25: Real part of the Hertz potential and of the field components of the localized 
electromagnetic ( "classic" , axially symmetric) X-shaped wave predicted, and first math- 
ematically constructed for the electromagnetic case, in refs.[T5]. For the meaning of 
the various panels, see the quoted references. The dimension of each panel is 4 m (in 
the radial direction) x 2 mm (in the propagation direction). [The values shown on the 
right-top corner of each panel represent the maxima and the minima of the images before 
normalization for display (MKSA units)] . 




Figure 26: Scheme of the experiment by Saari et al., who announced (PRL of 24 Nov. 1997) 
the production in optics of the beams depicted in the previous FigESl In the present figure 
one can see what it was shown by the experimental results: Namely, that the "X-shaped" 
waves are Superluminal: indeed, they, running after plane waves (the latter regularly 
travelling with speed c), do catch up with the considered plane waves. An analogous 
experiment has been later on performed with microwaves at Florence by Ranfagni et al. 
(PRL of 22 May 2000). 



Let us eventually recall the problem of producing an X-shaped Superluminal wave 
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Figure 27: In this figure a couple of elements are depicted of one of the trains of X-shaped 
pulses, mathematically constructed in ref . [113] , which propagate down a coaxial guide (in 
the TM case): This picture is just taken from ref. |113j . but analogous X-pulses exist (with 
infinite or finite total energy) for propagation along a cylindrical, normal-sized metallic 
waveguide] . 



like the one in Fig. 12, but truncated — of course- in space and in time (by use of a 
finite antenna, radiating for a finite time): in such a situation, the wave is known to 
keep its localization and Superluminality only till a certain depth of field [i.e., as long as 
they are fed by the waves arriving (with speed c) from the antenna], decaying abruptly 
afterwards. [36l [38] Let us add that various authors, taking account, e.g., of the time 
needed for fostering such Superluminal waves, have concluded that these localized Super- 
luminal pulses are unable to transmit information faster than c. Many of these questions 
have been discussed in what precedes; for further details, see the second of refs.|T5]. 

Anyway, the existence of the X-shaped Superluminal (or Super-sonic) pulses seem to 
constitute, together, e.g., with the Superluminality of evanescent waves, a confirmation of 
extended SR: a theory [51] based on the ordinary postulates of SR and that consequently 
does not appear to violate any of the fundamental principles of physics. It is curious 
moreover, that one of the first applications of such X-waves (that takes advantage of 
their propagation without deformation) has been accomplished in the field of medicine, 
and precisely — as we know — of ultrasound scanners [T9l [20] : while the most important 
applications of the (subluminal!) Frozen Waves will very probably affect, once more, 
human health problems like the cancerous ones. 

After the "digression" constituted by the above Appendix, let us go on to the Second 
Part of this work, with a slightly more technical [114] review about the physical and math- 
ematical characteristics of the Localized Waves and about some interesting applications. 
[In the Third Part we shall deal with the ones endowed with zero speed, i.e., with a static 
envelope, and, more in general, with the subluminal LWs]. 
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SECOND PART 



STRUCTURE OF THE NONDIFFRACTING WAVES 
AND SOME INTERESTING APPLICATIONS 



6 Foreword 

Since the early works [113- 116] on the so-called nondiffracting waves (or "Localized 
Waves"), a great deal of results has been published on this important subject, from 
both the theoretical and the experimental point of view. Initially, the theory was devel- 
oped taking into account only free space; however, in recent years, it has been extended 
for more complex media exhibiting effects such as dispersion [11 7- 11 9], nonlinearity [122] . 
anisotropy[l23j and losses p!2l] . Such extensions have been carried out along with the 
development of efficient methods for obtaining nondiffracting beams and pulses in the 
subluminal, luminal and Superluminal regimes [123- 130]. 

This Second Part addresses some theoretical methods related to nondiffracting solu- 
tions of the linear wave equation in unbounded homogeneous media, as well as to some 
interesting applications of such solutions. [1 14] 

The usual cylindrical coordinates (p, 0, z) will be used here. We already know that in 
these coordinates the linear wave equation is written as 



In Section VII we analyse the general structure of the Localized Waves, develop the so 
called Generalized Bidirectional Decomposition, and use it to obtain several luminal and 
Superluminal nondiffracting wave solutions of eq. fl27|) . 

In Section VIII we develop a kind of space-time focusing method by a continuous 
superposition of X-Shaped pulses of different velocities. 

Section IX addresses the properties of chirped optical X-Shaped pulses propagating 
in material media without boundaries. 

Subsequently, on the basis of what expounded in Section VII, we shall show at the 
beginning of the Third Part (in Section X) how a suitable superposition of Bessel beams 
can be even used to obtain stationary localized wave fields with high transverse local- 
ization, and whose longitudinal intensity pattern can assume any desired shape within a 
chosen interval < z < L of the propagation axis. 




(27) 
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For containing the length of this review, we had, obviously, to skip many interest- 
ing results. Let us just mention, for example, that rather simple analytic expressions, 
capable of describing the longitudinal (on-axis) evolution of axially-symmetric nondiffract- 
ing pulses, have been recently worked out in ref. [192j even for pulses truncated by finite 
apertures. By comparing what easily provided by such expressions, for several situa- 
tions (involving subluminal, luminal, or Superluminal localized pulses), with the results 
obtained by numerical evaluations of the Rayleigh-Sommerfeld diffraction integrals, an 
exellent agreement has been found. Therefore, those new closed-form expressions dis- 
pense with the need of time-consuming numerical simulations (and provide an effective 
tool for finding out the most important properties of the truncated localized pulses). 



7 Spectral structure of the Localized Waves and the 
Generalized Bidirectional Decomposition 

An effective way to understand the concept of the (ideal) nondiffr acting waves is furnishing 
a precise mathematical definition of these solutions, so to extract the necessary spectral 
structure from them. 

Intuitively, an ideal nondiffracting wave (beam or pulse) can be defined as a wave 
capable of maintaining indefinitely its spatial form (apart from local variations) while 
propagating. 

We can express such a characteristic by saying that a localized wave has to possess 
the property [124, 125] 

^(p,0,z,t) = ^(p,0,2; + Azo,t+^) (28) 

where Azq is a certain length and V is the pulse propagation speed that here can assume 
any value: < < oo. 

In terms of a Fourier Bessel expansion, we can write a function \l/(p, 0, z, t) as 



°° r roo roo roo ^ 

^'(p,0,z,t) = / dfcp / dk, / du kpA^{kp,k,,uj)Jnikpp)e'^^'t 



-iujt ^in<f> 



(29) 

On using the translation property of the Fourier transforms T[f{x+a)\ = exp{ika) T[f{x)], 
we have that A^{kp, k^, uj) and exp[i{kzAzo —ujAzo/V)]A^{kp, k^, uj) are the Fourier Bessel 
transforms of the l.h.s and r.h.s. functions in eq. (!28|) . And from this same equation we 
can get [124,125] the fundamental constraint linking the angular frequency u and the 
longitudinal wavenumber kz'- 



UJ 



Vkz + 2m7r 



V 



(30) 
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with m an integer. Obviously, this constraint can be satisfied by means of the spectral 
functions A^{kp, kz, u). 

Now, let us explicitly mention that constraint ( l30l) does not imply any breakdown of 
the wave equation. In fact, when inserting expression ( |29l) in the wave equation ( |271) . one 
gets that 



— - + 



(31) 



So, to obtain a solution of the wave equation by fl29|) . the spectrum A[^{kp, kz,uj) must 
possess the form 



A'^{kp,kz,uj) = An{kz,uj)6 



ki-r^-k 



(32) 



where is the Dirac delta function. With this, we can write a solution of the wave 
equation as 



uj/c 



duj / dkzAn{kz,uj)Jn p\ 



^2 I ^ik^z^-iut^ini 



(33) 



where we have considered positive angular frequencies only. 

Equation fl33l) is a superposition of Bessel beams and it is understood that the inte- 
grations in the (cu, /c^) plane are confined to the region < < oo and —u/c < kz <uj/c. 

Now, to obtain an ideal nondiffracting wave, the spectra An{kz.,uj) must obey the 
fundamental constraint (130|) . and so we write 



A„(/c^,cj) 



^ Snm{'^)5[uJ - {Vkz + K 



(34) 



where hm are constants representing the terms 2mT[Vl/S.ZQ in eq.( l30i) . and Snmiyj) are 
arbitrary frequency spectra. 

By inserting eq. llMI) into eq.( l33i) . we get a general integral form of the ideal non- 
diffracting wave fl28l) : 



^(p,0,2;,t) = Y Yl V^nm(p,0, 2;,t) 



(35) 



ra=— oo m=— oo 



with 



/■('^max)m 



(36) 
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where ujrnin and tUmax depend on the values of V: 

• for subluminal (V < c) locahzed waves: 6^ > 0, {ujjam)m = c6m/(c+ V^) and 

(t^max)m = cbm/{c - V)] 

• for luminal (V = c) localized waves: 6^ > 0, {oj^m)m = and {ujrna.x)m = cxo; 

• for Superluminal (V > c) localized waves: bm > 0, (co'min)m = c6m/(c + V^) and 

(t^max)m = OO. Or 6„ < 0, (u;min)m = cbm/{c- V) and (u;max)m = OO. 

It is important to notice that each ipnm{p, 4>j ^^t) the superposition ( 135|) is a truly 
nondiffr acting wave (beam or pulse) and the superposition of them, ( |35l) . is just the most 
general form to represent a nondiffracting wave defined by eq.( l28i) . Due to this fact, 
the search for methods capable of providing analytic solutions for ipnm{p, (p, z, t), eq. (|36i) . 
becomes an important task. 

Let us recall that equation (1361) is also a Bessel beam superposition, but with con- 
straint ( 130|) linking their angular frequencies and longitudinal wavenumbers. 

In spite of the fact that expression fl36l) represents ideal nondiffracting waves, it is 
difficult to obtain closed analytic solutions from it. Due to this, we are going to develop a 
method capable of overcoming such a difficulty, providing several interesting localized wave 
solutions (luminal and Superluminal) of arbitrary frequencies, including some solutions 
endowed with finite energy. 



7.1 The Generalized Bidirectional Decomposition 

For reasons that will be clear soon, instead of dealing with the integral expression ( l35l) . 
our starting point is the general expression ( |33l) . Here, for simplicity, we shall restrict 
ourselves to axially symmetric solutions, adopting the spectral functions 

An{kz,uj) = 6noMkz,uj) (37) 

where 6no is the Kronecker delta. 

In this way, we get the following general solution (considering positive angular fre- 
quencies only), which describes axially symmetric waves: 



^{p,(f),z,t) = du J_^^^dhA{h,u)Jo{p^-^-kl\e''''e-'^' (38) 

As we have seen, we can obtain ideal nondiffracting waves, given that the spectrum 
A{kz,u>) satisfies the linear relationship fl5U]) . Therefore, it becomes natural to choose 
new spectral parameters, in place of {uj, kz), that make easier to implement the mentioned 
constraint [124, 125]. With this in mind, let us choose the new spectral parameters (a,/3) 

a = ^{u + Vkz) ; f3 = ^(uj-Vk,). (39) 
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Let us consider here only luminal {V = c) and Superluminal {V > c) nondiffracting pulses. 

With the change of variables (|39l) in the integral solution ( l38i) . and considering {V > c), 
the integration limits on a and (3 have to satisfy the three inequalities 



< a + P < oo 



a > -P 



a > —p 



(40) 



c-V 

Let us suppose both a and j3 to be positive [a, /? > 0]. The first inequality in ( HOl) is then 
satisfied; while the coefficients {c — V)/{c+V) and {c + V)/{c — V) entering relations fHU]) 
are both negatives (since V > c). As a consequence, the other two inequalities in ( HOl) 
result to be automatically satisfied. In other words, the integration limits < a < oo 
and < /9 < oo are contained" into the limits (HOj) and are therefore acceptable. Indeed, 
they constitute a rather suitable choice for facilitating all the subsequent integrations. 

Therefore, instead of eq. fl38l) . we shall consider the (more easily integrable) Bessel 
beam superposition in the new variables [with > c] 



where we have defined 



oo roo 

da / d(3A(a, (3) Jo [ p ^ i 
Jo \ \ \ 



(a2 + /32) + 2 + 1 J ap \ e*°^e 



(41) 



C = z-Vt- r] = z + Vt . (42) 

The present procedure is a generahzation of the so-called "bidirectional decomposition" 
technique [125] . which was devised in the past for V = c. 

From the new spectral parameters defined in transformation (139|) . it is easy to see 
that the constraint ( !30l) . i.e. uj = Vkz + b, is implemented just by making 



A{h,u;)^A{a,P) = S{a)6{P - Po) (43) 

with Po = b/2V. The delta function 6{P — Po) in the spectrum ( H3|l means that we are 
integrating Bessel beams along the continuous line u = Vkz + 2l^/?o and, in this way, the 
function S{a) will give the frequency dependence of the spectrum: S{a) — >• S{uj/V — Po)- 

This method constitutes a simple, natural way for obtaining pulses with field concen- 
tration on p = and at ( = ^ z = Vt. 

Now, it is important to emphasize |127j that, when /3o > in ( H3ll . superposition ( l4T]l 
gets contributions from both backward and forward travelling Bessel beams, corresponding 
to the frequency intervals VPo < u < 2VPo (where kz < 0) and 2VPo < uj < oo (where 
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kz > 0), respectively. Nevertheless, we can obtain physical solutions when rendering 
the contribution of the backward-travelling components negligible, by choosing suitable 
weight functions S{a). 

It is also worth noticing that we adopted the new spectral parameters a and (3 just 
to obtain (closed-form) analytic localized wave solutions: The spectral characteristics 
of these new solutions can be brought into evidence by using transformations (l39l) and 
writing the corresponding spectrum in terms of the usual u and kz spectral parameters. 

In the following, we consider some cases with Po = and > 0. 



7.1.1 Closed analytic expressions describing some ideal nondiffracting pulses 



Let us first consider, in eq. fHT]) . spectra of the type fH3l) with Pq = 0: 

A{a,p) = aV6{p)e-^^" (44) 

A{a,p) = aV 5{p)Jo{2d^)e-~''^'' (45) 

A(a,/3) = 5{P)^^^e-'^^% (46) 
a 

a > and d being constants. 

One can obtain from the above spectra the following Superluminal LW solutions, 
respectively: 

— from spectrum (jSj), we can use the identity (6.611.1) in ref. |133] . obtaining the well 
known ordinary X xave solution (also called X-shaped pulse) 

*(P,C) ^ X = ""^^ , ' (47) 

^{aV - iCY + (f - l) 

— by using spectrum fH51) and the identity (6.6444) of ref. |133j . one gets 



/y2 



m{p,C) = X Jo\\l—-l {aVy'd'X'p ) exp -{aV - {aVy'd' X' ; (4J 



the Superluminal nondiffracting pulse 
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= sm 



-1 




X-2 + {d/aVf + 2pd{aV)-^^V^/c^ - 1 



1' 



(49) 



+ + {d/aVf - 2pd(aV)-y VVc2 - 1 



) 



is obtained from spectrum (146|) by using identity (6.752.1) of ref. |133 j for a > and d > 
From the previous discussion, we get to know that any solutions obtained from spectra 

of the type (l43l) with /5o = are free from noncausal (backward travelhng) components. 
In addition, when (3o = 0, we can see that the pulsed solutions depend on z and t 

through ( = z — Vt only, and so propagate rigidly, i.e. without distortion. Such pulses 

can be transversally localized only if > c, because if = c the function has to obey 

the Laplace equation on the transverse planes[124,125]. 

Many others Superluminal localized waves can be easily constructed [I27] from the 

above solutions just by taking the derivatives (of any order) with respect to (. It is also 

possible to show |127j that the new solutions, obtained in this way, have their spectra 

shifted towards higher frequencies. 

Now, let us pass to consider, in eq. pTl) . a spectrum of the type ( H3|) with /3o > 0: 



with a a positive constant. 

As we have seen, the presence of the delta function, with the constant (3q > 0, implies 
that we are integrating (summing) Bessel beams along the continuous line u = Vkz+2V (3q. 
Now, the function S{a) = aVexp{—aVuj) entails that we are considering a frequency 
spectrum of the type S{uj) oc exp(— acj), and therefore with a bandwidth given by Au = 
I /a. 

Since /?o > 0, the interval V[3q < oj < 2V (3q (or, equivalently, in this case, < a < /?o), 
corresponds to backward Bessel beams, i.e., negative values of kz- However, we can get 
physical solutions when making the contribution of this frequency interval negligible. In 
our case, this can be obtained by making a/3oV^ << 1, so that the exponential decay of the 
spectrum S with respect to u is very slow, and the contribution of the interval u > 2VPo 
(where kz > 0) largely overruns the VPo < u < 2VPo (where kz < 0) contribution. 

Incidentally, let us note that, once we ensure the causal behaviour of the pulse by 
making aVjSo << 1 in fl50|) . we have that Aa = 1/aV » Pq, and one can therefore 
simplify the argument of the Bessel function, in the integrand of superposition (HT]) . by 
neglecting the term (V'^/c^ — 1)i3q. With this, the superposition fHTj) . with the spectrum 
( 150|1 . can be written as 



A{a,P) 



aV6{(3 - Po) e 



—aVa 



(50) 
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C, v) ~ aV^e-^*'^ da Jo j p 



\ 



^ - 1 ) a' + 2 f ^ + 1 1 aPo 



(51) 



Now, we can use identity (6.616.1) of ref . |133j and obtain the new localized Superluminal 
solution called [T27] Superluminal Focus Wave Mode (SFWM): 



SFWM 



exp 



-C2 



{{aV - zO 



aVX- 



(52) 



where, as before, X is the ordinary X-pulse fHTI) . The center of the SFWM is located on 
p = and C = (i.e. at z = Vt). The intensity, |\l/|2, of this pulse propagates rigidly, it 
being a function of p and ( only. However, the complex function \E'sfwm (i-e. its real and 
imaginary parts) propagate with local variations, recovering their whole three dimensional 
form after each space and time interval Azq = 7r//3o and Ato = tt/(3oV. 

The SFWM solution written above, for V — > reduces to the well known Focus 
Wave Mode (FWM) solution [125], travelling with speed c: 



FWM 



ac- 



ac 



exp 



Pop' 



ac 



(53) 



Let us also emphasize that, since Pq > 0, spectrum (150|) results to correspond to angular 
frequencies lo > VPo. Thus, our new solution can be used to construct also high frequency 
pulses. 



7.1.2 Finite energy nondiffracting pulses 



In this subsection, we shall show how to get finite energy localized wave pulses, that can 
propagate for long distances while maintaining their spatial resolution, i.e., that possess 
a large depth of field. 

As we have seen, ideal nondiffracting waves can be constructed by superposing Bessel 
beams [cf. eg. (1381) for cylindrical symmetry] with a spectrum A{uj,kz) that satisfies a 
linear relationship between u) and kz. In the general bidirectional decomposition method, 
this can be obtained by using spectra of the type (H3|l in superposition (14T]) . 

Solutions of that type possess an infinity depth of field: however, they are endowed 
with infinite energy [123, 125]. To overcome this problem, we can truncate an ideal non- 
diffracting wave by a finite aperture, and the resulting pulse will have finite energy and 
a finite field-depth. Even so, such field-depths may be very large when compared with 
those of ordinary waves. 
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The problem in the present case is that the resulting field has to be calculated by 
diffraction integrals (such as the well known Rayleigh-Sommerfeld formula) and, in gen- 
eral, a closed analytic formula for the resulting pulse cannot be obtained. 

However, there is another way to construct localized pulses with finite energy [127]. 
Namely, by using spectra A{uj, k^) in eq. fISS]) whose domains are not restricted exactly to 
the straight line u = Vkz + b, but are defined in the surroundings of that line, wherein the 
spectra should have their main values concentrated (in other words, any spectrum has to 
be well localized in the vicinity of that line). 

Similarly, in terms of the generalized bidirectional decomposition given in eq. pTl) . 
finite energy nondiffracting wave pulses can be constructed by adopting spectral functions 
A{a,P) well localized in the vicinity of the line P = Pq, quantity Pq being a constant. 

To exemplify this method, let us consider the following spectrum 



aqV e' 



-aVa^-q{l3-f3o) 



for P>Po 



A{a,P) 



(54) 



for < /? < /?o 

in superposition fHTl) . quantities a and q being free positive constants and V the peak's 
pulse velocity (here V > c). 

It is easy to see that the above spectrum is zero in the region above the P = Pq line, 
while it decays in the region below (as well as along) such a line. We can concentrate this 
spectrum on P = Pq by choosing values of q in such a way that qPo >> 1. The faster the 
spectrum decay takes place in the region below the P = Pq line, the larger the field-depth 
of the corresponding pulse results to be. 

Besides this, once we choose qPo >> 1 to obtain pulses with a large field-depth, we 
can also minimize the contribution of the noncausal (backward) components by choosing 
aVPo << 1; analogously with what we did and obtained for the SFWM case. 

Again in analogy with the SFWM case, when we choose qPo » 1 (i.e., a long field- 
depth) and aVPo << 1 (a minimal contribution of the backward components), one is 
allowed to simplify the argument of the Bessel function, in the integrand of superposition 
( liTl) . by neglecting the term [V'^/c^ — l)Po. 

With the help of the observations above, one can write the superposition ( HTi) . with 
spectrum (IMj), as: 



roo roo I /y 

*(p,C,r/) ^ agl^y^ dp daJolp^)^^ 



— -l] + 2(^ + 1] ap 



and, using identity (6.616.1) given in ref. |133] . we get 



roo 

^{p,C,v) -qx / dpe 



-(j(/3-/3o)g-j/3»?gjjp 



p 



y2 _c2 



(aV -iC- aVX~^ 



(55) 



(56) 
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which can be viewed as a superposition of the SFWM pulses (see eq. flS21) ). 

The above integration can be easily performed and results [T2Zj in the so called Su- 
perluminal Modified Power Spectrum (SMPS) pulse: 



,T, ( /- \ ^ exp[(y-i7/)/?o] 

^sMPs(p,C,^) = — — — - 




(57) 



where X is the ordinary X-pulse P7|) and Y is defined by 



Y 



y2_c2 - 



{aV - i() - aVX 



-1 



(58) 



The SMPS pulse is a Superluminal localized wave, with field concentration around p = 
and C = (i-6-, z = Vt), and with finite total energy. We will show that the depth of 
field, Z, of this pulse is given by Z^mps = Q'/2. 

An interesting property of the SMPS pulse is related to its transverse width (the 
transverse spot size at the pulse center). It can be shown from eq.( l57|) that, for the cases 
where aV « l/Po and g/3o >> 1, i.e., for the cases considered by us, the transverse spot 
size, Ap, of the pulse center {( = 0) is determined by the exponential function in (1571) 
and is given by 



which clearly does not depend on z, and so remains constant during its propagation. In 
other words, in spite of the fact that the SMPS pulse suffers an intensity decrease during 
propagation, it preserves its transverse spot size. This interesting characteristic is not 
met in ordinary pulses, like the gaussian ones, where the amplitude of the pulse decreases 
and the width increases by the same factor. 

Figure [28] shows the intensity of a SMPS pulse, with Pq = 33 m^^, V = 1.01c, a = 
10~^^ s and q = 10^ m, at two different moments, for t = and after 50 km of propagation, 
where, as one can see, the pulse becomes less intense (precisely, with half of its initial 
peak intensity). In spite of the intensity decrease, the pulse maintains its transverse 
width, as one can see from the 2D plots in Fig. (!28ll . which show the field intensities in 
the transverse sections at 2; = and z = q/2 = 50 km. 

Other three important well known finite energy nondiffr acting solutions can be ob- 
tained directly from the SMPS pulse: 

— The first one, obtained from ( 1571) by making /?o = 0, is the so called [T27,j Superluminal 
Splash Pulse (SSP) 



Ap = c 



aV - c2 



(59) 



2 ' 



^ssp(p,C,^) 



qX 



(60) 



q + irj — Y 
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Figure 28: Representation of a Superluminal Modified Power Spectrum pulse, eg. ( 1571) . 
Its total energy is finite (even without any truncation), and so it gets deformed while 
propagating, since its amplitude decreases with time. In Figj28h we represent, for t = 0, 
the pulse corresponding to /3o = 33 m~^, V = 1.01c, a = 10~^^s and q = 10^ m. In 
FigJ25b it is depicted the same pulse after having travelled 50 km. 



— The other two pulses are luminal. By taking the limit — c"*" in the SMPS pulse 
(|57|) . we get the well known [125] luminal Modified Power Spectrum (MPS) pulse 

^"^^^^'^'^^ = (g + zr,)(ac-zC)+p^ ^"P K^C) ' ^''^ 

finally, by taking the limit V and making /3o = in the SMPS pulse [or, equivalently, 

by making /3o = in the MPS pulse (|6T|) . or, instead, by taking the limit \^ ^ c"*" in the 
SSP (l60l) ]. we obtain the well known[T25] luminal Splash Pulse (SP) solution 

^sp(p,c,^) = T—-zr^^-7v—^- (62) 

[q + iri)[ac - iQ + 

It is also interesting to notice that the X and SFWM pulses can be obtained from the SSP 
and SMPS pulses (respectively) by making g — > oo in Eqs. (1601) and ( 1571) . As a matter of 
fact, the solutions SSP and SMPS can be viewed as the finite energy versions of the X 
and SFWM pulses, respectively. 



Some characteristics of the SMPS pulse: 

Let us examine the on-axis (p = 0) behaviour of the SMPS pulse. On p = we have 

*SMPs(p = 0,C,r7) = aqVe-'C^'liaV -iC){q + ir,)]-^ . (63) 

From this expression, we can show that the longitudinal localization A^, for t = 0, of the 
SMPS pulse square-magnitude is 
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= 2aV . (64) 

If we now define the field-depth Z as the distance over which the pulse's peak intensity 
is at least 50% of its initial valu^, then we can obtain, from (jHS]), the depth of field 

^SMPS = I ! (65) 

which depends only on g, as we expected since q regulates the concentration of the spec- 
trum around the line uj = Vk^ + ^VPo. 

Now, let us examine the maximum amplitude M of the real part of fl551) . which for 
z = Vt writes (C = and r] = 2z): 

MsMPs = Re WsMPs(p = 0,z = Vt)\ = , , , ^^ • (66) 

1 + ^{z/qY 

Initially, for z = 0, t = 0, one has M = 1 and can also infer that: 

(i) when z/q << 1, namely when z << Z, equation (!66l) becomes 

MsMPs ~ cos(2/?o^) for z « Z (67) 

and the pulse's peak actually oscillates harmonically with "wavelength" Azq = 7i/j3o and 
"period" Ato = 'n'/VPo, all along its field-depth; 

(ii) when z/q » 1, namely z » Z, equation ([HS]) becomes 

MsMPS « -'^^^^^ for z»Z. (68) 
2z/q 

Therefore, beyond its depth of field, the pulse goes on oscillating with the same Azq, but 
its maximum amplitude decays proportionally to z. 

In the next two Sections we are going to see applications of the localized wave pulses. 



8 Space-Time Focusing of X-shaped Pulses 

In this Section we are going to show how one can in general use any known Superluminal 
solution to obtain from it a large number of analytic expressions for space-time focused 
waves, endowed with a very strong intensity peak at the desired location. The method 
presented here is a natural extension of that developed by A.Shaarawi et al. [134] . where 

^We can expect that, while the pulse peak-intensity is maintained, the pulse keeps its spatial form 
too. 
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the space-time focusing was achieved by superimposing a discrete number of ordinary 
X-waves, characterized by different values 6 of the axicon angle. 

In this Section, based on ref. [T35], we shall go on to more efficient superpositions for 
varying velocities V, related to 9 through the known[115, 116,54] relation V = c/ cos6. 
This enhanced focusing scheme has the advantage of yielding analytic (closed-form) ex- 
pressions for the spatio-temporally focused pulses. 

Let us start by considering an axially symmetric ideal nondiffracting Superluminal 
wave pulse ip{p, z — Vt) in a dispersionless medium, where V = c/ cos 9 > c is the pulse 
velocity, 9 being the axicon angle. As we have seen in the previous Section, pulses like 
these can be obtained by a suitable frequency superposition of Bessel beams. Suppose 
that we have now N waves of the type z — Vn{t — tn)), with different velocities 

c < Vi < V2 < .. < Vjy, and emitted at (different) times t„; quantities t„ being constants, 
while n = 1,2, ...N. The center of each pulse is located ai z = Vn{t — t„). To obtain a 
highly focused wave, we need all the wave components ipn{p, z — Vnit — tn)) to reach the 
given point, z = Z{, at the same time t = tf. On choosing ti = for the slowest pulse ipi, 
it is easily seen that the peak of this pulse reaches the point z = Zf at the time tf = Z{/Vi. 
So we obtain that, for each ipn, the instant of emission tn must be 

*„ = (i - ^) . (69) 

With this in mind, we can construct other exact solutions to the wave equation, given by 

= C '^''■^c)^''-'- - " (' - (v^ 4) ^')) ' (™) 

where V is the velocity of the wave ip{p.,z — Vt) which enters the integrand of (ITOll . 
While integrating, V is considered as a continuous variable in the interval [Vmin, ^ax]- In 
eq.(70), function A{y) is the velocity distribution that specifies the contribution of each 
wave component (with velocity V) to the integration. The resulting wave \E'(p, ^,t) can 
have a more or less strong amplitude peak aX z = Z{, at time tf = Zf/Vmm, depending on 
A{V) and on the difference Knax — Knin- Let us notice that also the resulting wavefield will 
propagate with a Superluminal peak velocity, depending on A(y) too. In the cases when 
the velocity-distribution function is well concentrated around a certain velocity value, 
one can expect the wave (|70i) to increase its magnitude and spatial localization while 
propagating. Finally, the pulse peak acquires its maximum amplitude and localization 
at the chosen point z = Zf, and at time t = Zf/Vj^m, as we know. Afterwards, the wave 
suffers a progressing spreading, and a decreasing of its amplitude. 

8.1 Focusing Effects by Using Ordinary X- Waves 

Here, we present a specific example by integrating ( ITOl) over the standard, classic pT8] 
X-waves, X = aV[{aV — i{z — Vt))'^ + {V'^/c^ — l)p^]~^/^. When using this ordinary 
X-wave, the largest spectral amplitudes are obtained for low frequencies. For this reason. 
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one may expect that the solutions considered below will be suitable mainly for low fre- 
quency applications. Let us choose, then, the function ip in the integrand of eq. (J70l) to be 
^(p, t) = X{p, z-V{t- (1/Kni„ - l/V)z,)), viz.: 



aV 



^[aV-,(z-V(t- - i) .,))r + (f - 1) P'- 

After some manipulations, one obtains the analytic integral solution 



^(p,z,t) 



aVA{V) 



y^i„ ^yPV + QV + R 



dV 



with 



(71) 



(72) 



P 



a + i It 



^^min 



+ 77 



Q = 2(t- 



Vmin 



ai^ {z — Zf) 



R = [-{z - Zi) 



(73) 



In what follows, we illustrate the behaviour of some new spatio-temporally focused 
pulses, by taking into consideration a few different velocity distributions A{V). These 
new pulses are closed analytic exact solutions of the wave equation. 



First example: 

Let us consider our integral solution fl72l) with A{V) = 1 s/m. In this case, the contri- 
bution of the X-waves is the same for all velocities in the allowed range [\4iin, Knax]- On 
using identity 2.264.2 listed in ref. |133j . we get the particular solution 

^(P,^,i) = I (V^K^ax + QVm.. + R- ^PV^^ + QV^^ + R) 

aQ / 2 ^PjPVl,^ + QVmin + R)+ 2PKnin + Q \ ^^^^^ 

where P, Q and R are given in eq.( 173i) . A 3-dimensional (3D) plot of this function is 
provided in Figl29l where we have chosen a = 10~^^ s, Vmin = 1.001 c, \4iax = 1.005 c 
and Zf = 200 cm. It can be seen that this solution exhibits a rather evident space-time 
focusing. An initially spread-out pulse (shown for t = 0) becomes highly localized at 
t = t{ = 2;f/Ktiin = 6.66 ns, the pulse peak amplitude at Z{ being 40.82 times greater than 
the initial one. In addition, at the focusing time tf the field is much more localized than 
at any other times. The velocity of this pulse is approximately V = 1.003 c. 



54 




Figure 29: Space-time evolution of the Superluminal pulse represented by eq-f TM]) : the 
chosen parameter values are a = 10~^^ s; V^nin = 1-001 c; Vmax = 1-005 c while the 
focusing point is at Zf = 200 cm. One can see that this solution is associated with 
a rather good spatio-temporal focusing. The field amplitude aX z = Zf is 40.82 times 
larger than the initial one. The field amplitude is normalized at the space-time point 
p = 0, z = Zf, t = tf. 



Second example: 

In this case we choose A{V) = 1/V s/m, and, using the identity 2.261 in ref. |133j . eq. (!72|) 
yields 

a , (2 ^P{PVl,^ + QKnax + R)+ 2PK.ax + Q\ 

^{p,z,t) = In V . (75) 

VP \ 2 ^P{PV^^ + gKnin + R)+ 2PKni, + Q I 

Other exact closed-form solutions can be obtained jl35j by considering, for instance, ve- 
locity distributions like A{V) = and A{V) = 

Once more, we can actually construct many others spatio-temporally focused pulses 
from the above solutions, just by taking their time derivatives (of any order). It is possible 
to show [135] that also the new solutions obtained in this way have their spectra shifted 
towards higher frequencies. 
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9 Chirped Optical X-Type Pulses in Material Media 



The theory of the locahzed waves was initially developed for free space (vacuum). In 
1996, SSnajalg et al. |119j showed that the localized wave theory can be extended to include 
(unbounded) dispersive media. This was obtained by making the axicon angle of the Bessel 
beams (BBs) vary with the frequency[117-119] in such a way that a suitable frequency 
superposition of these beams does compensate for the material dispersion. Soon after this 
idea was reported, many interesting nondiffracting/nondispersive pulses were obtained 
theoretically[117-119] and experimentally [1 19j . In spite of such an extended method to 
be of remarkable importance, working well in theory, its experimental implementation is 
not so simpl^. 

In 2004 Zamboni-Rached et al. [136] developed a simpler way to obtain pulses capable 
of recovering their spatial shape, both transversally and longitudinally, after some prop- 
agation. It consisted in using chirped optical X-typed pulses, while keeping the axicon 
angle fixed. Let us recall that, by contrast, chirped Gaussian pulses in unbounded ma- 
terial media may recover only their longitudinal shape, since they undergo a progressing 
transverse spreading while propagating. 

The present Section is devoted to this approach. Let us start with an axially- 
symmetric Bessel beam in a material medium with refractive index n{uj): 

ilj{p,z,t) = jQ{kpp) exp{iPz) exp{—iLjt) , (76) 

where it must be obeyed the condition = r?[LS)iJ^ j (? — which connects among 
themselves the transverse and longitudinal wave numbers kp and /5, and the angular 
frequency lo. In addition, we impose that A;^ > and > 0, to avoid a nonphysical 
behaviour of the Bessel function Jo(-) cind to confine ourselves to forward propagation only. 
Once the conditions above are satisfied, we have the liberty of writing the longitudinal 
wave number as /? = {n[LS)uj co-^ff) j c and, therefore, kp = {n{uj)uj sin6)/c; where (as in 
the free space case) 6 is the axicon angle of the Bessel beam. 

Now we can obtain an X-shaped pulse by performing a frequency superposition of 
these Bessel beams [BB], with /3 and kp given by the previous relations: 

S{uj) Jo I ^ ^ singpj exp[i(3{uj)z] exp{-iujt) du , (77) 

where S{u) is the frequency spectrum, and the axicon angle is kept constant. One can 
see that the phase velocity of each BB in our superposition fl77p is different, and given by 
Vphase = c/{n{ijj) COS 6). So, the pulse represented by eq. flTTl) will suffer dispersion during 
its propagation. 

llWe refer the interested reader to quotations[117-119] for obtaining a description, theoretical and 
experimental, of that extended method 
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As we said, the method developed by Sonajalg et al. [119] . and explored by oth- 
ers [11 8, 11 9], to overcome this problem consisted in regarding the axicon angle 6* as a 
function of the frequency, in order to obtain a linear relationship between f3 and u. 

Here, however, we choose to work with a fixed axicon angle, and we have to find 
out another way for avoiding dispersion and diffraction all along a certain propagation 
distance. To do that, we might choose a chirped gaussian spectrum ^(a;) in eq. (|771) 



S{u;) 



Tn 



exp[-g [uj-ujo) 



with 



'T'2 
^0 



2{1 + iC) 



(78) 



/27r(l 

where cuq is the central frequency of the spectrum, Tq is a constant related with the 
initial temporal width, and C is the chirp parameter (we chose as temporal width the 
half-width of the relevant gaussian curve when its heigth equals 1/e times its full heigth). 
Unfortunately, there is no analytic solution to eq.( 177|) with S{uj) given by eq.llTHl). so that 
some approximations are to be made. Then, let us assume that the spectrum ^(ci;), 
in the surroundings of the carrier frequency Uq , is narrow enough to guarantee that 
Auj/uq << 1, so to ensure that l3{uj) can be approximated by the first three terms of 
its Taylor expansion in the vicinity of Uq: That is, P{uj) ^ P{uJo) + P'{u!)\^q {u — ujq) + 
(l/2)/5"(u;)|tj„ (cu — uqY] when, after using /5 = n{Lj)uj cos 9 / c, it results that 



duj 



cos 6 



. , dn 
n{uj) +UJ — 



cos 9 



2 dn d^n 



(79) 



As we know, 13' {u) is related to the pulse group- velocity by the relation Vg = 
Here we can see the difference between the group- velocity of the X-type pulse (with a fixed 
axicon angle) and that of a standard gaussian pulse. Such a difference is due to the factor 
cos 9 in eq. fl79|) . Because of it, the group- velocity of our X-type pulse is always greater than 



the gaussian's. In other words, (V^)x = (V cos6')(Vg)gauss- We also know that the second 
derivative of P{uj) is related to the group- velocity dispersion (GVD) P2 by P2 = 
The GVD is responsible for the temporal (longitudinal) spreading of the pulse. Here one 
can see that the GVD of the X-type pulse is always smaller than that of the standard 
gaussian pulses, due the factor cos6' in eq. (j79l) . Namely: {P2)x = cos 6'(/32) gauss- 
On using the above results, we can write 



Tq exp[iP{uJo)z] exp{—iuJot) 

pTx{l + iC) 



duj Jo 



n[Uii)uj 



sin 9 p 



(80) 



X exp < z- 



.[U - UIq) 



[z — Vgt] > exp — ouqY 



iP2 2 



The integral in eq. flHOl) cannot be evaluated analytically, but for us it is sufficient to obtain 
the pulse behaviour. Let us analyse the pulse for p = 0. In this case we get 
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^(p = 0,z,t) 



To exp[iP{uJo)z] exp{—iuot) 



exp 



-{z-Vgtf{l + iC) 



(81) 



2V^^[Ti-if32il + tC)z] 



From eq. flHTl) one can immediately see that the initial temporal width of the pulse intensity 
is To and that, after a propagation distance z, the time-width Ti becomes 



Relation fl82l) describes the pulse spreading-behaviour. One can easily show that such a 
behaviour depends on the sign (positive or negative) of the product l32C, as is well known 
to happen for the standard gaussian pulses ^137] . In the case P2C > 0, the pulse will 
monotonically become broader and broader with the distance z. On the other hand, if 
P2C < the pulse will suffer, in a first stage, a narrowing, and then (during the rest of 
its propagation) it will spread . So, there will be a certain propagation distance AT after 
which the pulse will recover its initial temporal width (Ti = To). From relation ( l82l) . we 
can find such a distance Zti=To (considering < 0) to be 



One may notice that the maximum distance at which our chirped pulse, with given Tq 
and /?2, may recover its initial temporal width can be easily evaluated from eq. (j83l) . and 
it results to be Tdisp = T^/P2- We shall call such a maximum value Tdisp the "dispersion 
length": It is the maximum distance the X-type pulse may travel while recovering its 
initial longitudinal shape. Obviously, if we want the pulse to reassume its longitudinal 
shape at some desired distance z < L^isp, we have just to suitably choose a new value for 
the chirp parameter. 

Let us emphasize that the property of recovering its own initial temporal (or longitu- 
dinal) width may be verified to exist also in the case of chirped standard gaussian pulses. 
However, the latter will suffer a progressing transverse spreading, which will not be re- 
versible. The distance at which a gaussian pulse doubles its initial transverse width wq 
is Zdis = V^t^Wq/ Xq, where Ao is the carrier wavelength. Thus, one can see that optical 
gaussian pulses with great transverse localization will get spoiled in a few centimeters or 
even less. 

Now we shall show that it is possible to recover also the transverse shape of the 
chirped X-type pulse intensity; actually, it is possible to recover its entire spatial shape 
after a distance Zti=To- To see this, let us go back to our integral solution (IHOl) . and 
perform the change of coordinates {z,t) {Az,tc = Zc/Vg), with 




(82) 



(83) 
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Z = Zc + IS.Z 



t = U = 



where Zc is the center of the pulse {Az is the distance from such a point), and tc is the 
time at which the pulse center is located at Zc- What we are going to do is comparing 
our integral solution (IHOj) . when Zc = (initial pulse), with that when Zc = -^Ti=ro = 
—2CT^/{[32{C-2 + 1)). In this way, solution flHOj) can be written, when Zc = 0, as 



^{p,z, = 0,Az) 



Tq exp{iPoAz) 



did jQ{kp{ijj)p) exp 



2(1 + ^2) 



^5) 



X exp < I 



{to - ujq)Az_ ^ {uj - uJoyi32Az ^ (cu - ujoYT^C 



2(1 + 



where we have taken the value q given by eq. 
recovers its entire original form at Zc = Zti=To 
our integral solution at that point, obtaining 



To verify that the pulse intensity 
-2CToV[/?2(C2 + 1)], we can analyse 



To exp <^ i/3o 



^(p, = ^Ti=To, Az) 



Zr - Az' 



CZr, 



cos6n{(jJo)Vg 



2n{l + iC) 



X / duj Jo{kp{uj)p) exp 



2(1 + C2) 



(86) 



X exp < —t 



2 2(1 + C2) 



where we put Az = —Az'. In this way, one immediately sees that 



\'^ip,z, = 0,Az)\^ = \^ip,Zc = Zt,=To,-^z) 



17) 



59 



Therefore, from eq. fl57|) it is clear that the intensity of a chirped optical X-type pulse 
is able to recover its original three-dimensional shape, just with a longitudinal inversion 
at the pulse center. The present method results to be, therefore, a simple and effective 
procedure for compensating diffraction and dispersion in an unbounded material medium; 
and a method simpler than the one of varying the axicon angle with the frequency. 

Let us stress again that one can determine the distance z = ^Ti=To < -Z^disp at which 
the pulse takes on again its spatial shape by choosing a suitable value of the chirp param- 
eter. 

We have shown that the chirped X-type pulse recovers its three-dimensional shape 
after some distance, and we have also obtained an analytic description of the pulse longi- 
tudinal hehayioui (for p = 0) during propagation, by means of eq. llHT]) . However, we have 
not got yet the same information about the pulse transverse behaviour: We just learned, 
till now, that it will be recovered at z = Zt^=To- 

So, to complete the picture, we should find out also the transverse behaviour in the 
plane of the pulse center z = Vgt: We would then obtain quantitative information about 
the evolution of the pulse-shape during its entire propagation. But we are not going to 
expound all the relevant mathematical details here; let us only state that the transverse 
behaviour of the pulse (in the plane z = Zc = Vgt), during its whole propagation, can 
approximately be described by 



^ip,z 



Zci t 



Zc/Vg, 



To exp[i/9(co'o)-2] exp(— icjot) 
2ti{1 + iC) 



exp 



- tan^ e 



X 



r(i/2)Jo 



' n{ijjQ) ujq sin 6 



tan^ e p2 



p=l 



r(2j9 + i) 



'2p 



where /p(.) is the modified Bessel function of the first kind of order p, quantity r(.) being 



the gamma function and q being given by flTSl) . The interested reader can check ref. |136j 
for details on how eq. (l88!) is obtained from eq. flSOl) . 

At a first sight, this solution appears to be very complicated, but the series in its r.h.s. 
gives a negligible contribution. This circumstance renders our solution (!88|) of important 
practical interest, and we will use it in the following. For additional information about 
the transverse pulse evolution (to be extracted from eq. (l88|l ). the reader can consult again 
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ref . |136] . In the same paper, it is analysed how the generation by a finite aperture affects 
the chirped X-type pulses. 

The valuable methods developed in ref. [136j . and that we are partially revisiting in 
this Section, are of general interest, and work is in progress for applying them, e.g., also 
to the (different) case of the Schroedinger equation. 



9.1 An example: Chirped Optical X- typed pulse in bulk fused 
Silica 



For a bulk fused Silica, the refractive index n{uj) can be approximated by the Sellmeier 
equation |137] 



N R , ,2 

n\.) = 1 + E , (89) 

i=i '^i - 

where Uj are the resonance frequencies, Bj the strength of the j-th resonance, and N the 
total number of the material resonances that appear in the frequency range of interest. 
For our purposes it is appropriate to choose = 3, which yields, for bulk fused silica |137] . 
the values Bi = 0.6961663; B2 = 0.4079426; ^3 = 0.8974794; Ai = 0.0684043 /i m; 
A2 = 0.1162414 /i m; and A3 = 9.896161 m. 

Now, let us consider in this medium a chirped X-type pulse, with Aq = 0.2 /xm. 
To = 0.4 ps, C = —1, and with axicon angle 6 = 0.00084 rad: That correspond to 
an initial central spot with Apo = 0.117 mm. We get from eqs. flHTj) and fl88|) the 
longitudinal and transverse pulse evolution, which are represented in Fig|30l 

From Fig.(!30]a). we can observe that the pulse suffers initially a longitudinal narrowing 
with an increase of intensity, till the position z = Tq /2/32 = 0.186 m. After that point, the 
pulse starts broadening and decreasing its intensity, while recovering its entire longitudinal 
shape (width and intensity) at the point z = Tq / P2 = 0.373m, as it was predicted. At 
the same time, from Fig. (1301 b). one can notice that the pulse maintains its transverse 
width Ap = 2.4c/(n(a;o)a;osin6') = 0.117 mm (because Tqljo >> 1) during its entire 
propagation. The same does not occur, however, with the pulse intensity: Initially, the 
pulse suffers an increase of intensity, till position Zc = Tq/2/52 = 0.186 m; after that point 
the intensity starts decreasing, and the pulse recovers its entire transverse shape at point 
Zc = Tq / P2 = 0.373m, as expected. In the calculations we could skip the series in the 
r.h.s. of eq. (!88|) . because, as we already said, it yields a negligible contribution. 

Summarizing, from Figl30l we can see that the chirped X-type pulse recovers totally 
its longitudinal and transverse shape at position z = Ldisp = Tq / P2 = 0.373 m, as we 
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Figure 30: (a): Longitudinal-shape evolution of a chirped X-type pulse, propagating in 
fused sihca with Aq = 0.2/i m, Tq = 0.4 ps, C = —1, and axicon angle 9 = 0.00084 rad, 
which correspond to an initial transverse width of Apo = 0.117 mm. (b): Transverse- 
shape evolution for the same pulse. 

expected. Let us recall that a chirped gaussian pulse may just recover its longitudinal 
width, but with an intensity decrease, at the position given hy z = ^Ti=ro = -^^disp = 
Tq/132- Its transverse width, on the other hand, suffers a progressing and irreversible 
spreading. 

In the following Third Part we are going to "complete" our review by investigating 
also the (not less interesting) case of the suhluminal Localized Solutions to the wave 
equations, which, among the others, will allow us to set forth remarkable considerations 
about the role of (extended) Special Relativity. For instance, the various Superluminal 
and subluminal LWs are expected to be transformed one into the other by suitable Lorentz 
transformations. We shall start by studying, in terms of various different approaches, the 
very peculiar topic of zero-speed waves: Namely, the question of constructing localized 
fields with a static envelope; consisting, for example, in "light at rest" endowed with zero 
peak- velocity. We called Frozen Waves such solutions: They can have a lot of applications. 
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THIRD PART 



"FROZEN WAVES", 
AND THE SUBLUMINAL WAVE-BULLETS 



10 Modeling the Shape of Stationary Wave Fields: 
Frozen Waves 

As just mentioned, we start this Third Part by studying the very pecuhar topic of zero- 
speed waves: Namely, the question of constructing locahzed fields with a static envelope 
(for example, consisting in "light at rest" endowed with null peak- velocity). We called 
Frozen Waves such solutions: They permit a priori a lot of applications, as we are going 
to see. 

In the present Section we develop a very simple first method[122,128,129], based on 
our Second Part, by having recourse to superpositions of forward propagating and equal- 
frequency Bessel beams, that allows one controlling the longitudinal beam-intensity shape 
within a chosen interval < z < L, where z is the propagation axis and L can be much 
greater than the wavelength A of the monochromatic light (or sound) which is being 
used. Inside such a space interval, indeed, we succeed in constructing a stationary enve- 
lope whose longitudinal intensity pattern can approximately assume any desired shape, 
including, for instance, one or more high-intensity peaks (with distances between them 
much larger than A); and which — in addition — results to be naturally endowed also with 
a good transverse localization. Since the intensity envelopes remains static, i.e., with 
velocity V = 0, we called "Frozen Waves" (FW) such new solutions[122,128,129] to the 
wave equations. 

Although we are dealing here with exact solutions of the scalar wave equation, vec- 
torial solutions of the same kind for the electromagnetic field can be worked out: Indeed, 
solutions to Maxwell's equations may be naturally inferred even from the scalar wave- 
equation solutions[136-138]. 

We present first the method referring to lossless media[128,129] while, in the second 
part of this Section, we extend the method to absorbing media [T2^ . 
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10.1 Stationary wavefields with arbitrary longitudinal shape in 
lossless media, obtained by superposing equal-frequency 
Bessel beams 

Let us start from the well-known axis-symmetric zeroth order Bessel beam solution to the 
wave equation: 

ij{p,z,t) = JolMe'^^'e""* (90) 

with 

^2 = ^-/5^ (91) 

where u, kp and (3 are the angular frequency, the transverse and the longitudinal wave 
numbers, respectively. We also impose the conditions 

Lu/P > and > (92) 

(which imply uj/f3 > c) to ensure forward propagation only (with no evanescent waves), 
as well as a physical behaviour of the Bessel function Jq- 

Now, let us make a superposition of 2N + 1 Bessel beams with the same frequency 
ujq, but with different (and still unknown) longitudinal wave numbers jSm'- 

N 

^{p,z,t) = e-*-«* AmMkp^p)e'^-\ (93) 

rn=-N 

where the m represent integer numbers and the are constant coefficients. For each 
m, the parameters uq, kpm and Pm must satisfy eq. (l9Tl) . and, because of conditions (l92l) . 
when considering c^o > we must have 

0<Pm< — . (94) 
c 

Let us now suppose that we wish |\l/(p, t)p, given by eq. (p3|) . to assume on the axis 
p = the pattern represented by a function |F(z)p, inside the chosen interval < z < L. 
In this case, the function F{z) can be expanded, as usual, in a Fourier series: 



oo 

F{z) = Y Bme'"^ 

m=—oo 

where 



LJo 
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More precisely, our goal now is finding out the values of the longitudinal wave numbers 
and the coefficients Am of (jHS]), in order to reproduce approximately, within the said 
interval < z < L (for p = 0), the predetermined longitudinal intensity-pattern 
Namely, we wish to have 



N 
m=-N 



\F{z)\^ with 0<z<L. (95) 



Looking at eg. (1951) . one might be tempted to take Pm = 27im/L, thus obtaining a 
truncated Fourier series, expected to represent approximately the desired pattern F{z). 
Superpositions of Bessel beams with (3^, = 27rm/L have been actually used in some works 
to obtain a large set of transverse amplitude profiles p30] . However, for our purposes, this 
choice is not appropriate, due to two principal reasons: 1) It yields negative values for 
Prn (when m < 0), which implies backward propagating components (since ujq > 0); 2) 
In the cases when L » Aq, which are of our interest here, the main terms of the series 
correspond to very small values of (3rn, which results in a very short field-depth of the 
corresponding Bessel beams (when generated by finite apertures), preventing the creation 
of the desired envelopes far form the source. 

Therefore, we need to make a better choice for the values of (3rn, which permits forward 
propagation components only, and a good depth of field. This problem can be solved by 
putting 

Pm = Q+'^m, (96) 

where Q > is a value to be chosen (as we shall see) according to the given experimental 
situation and the desired degree of transverse field localization. Due to eq. (lM|) . one gets 

0<Q±—N<'^. (97) 
L c 

Inequality (1971) . can be used to determine the maximum value of m, that we call A^max, 
once Q, L and Uq have been chosen. 

As a consequence, for getting a longitudinal intensity pattern approximately equal to 
the desired one, |F(^)p, in the interval < z < L, eq. fl93|) has to be rewritten as 



N 

^^f{p = 0,z,t) = e-'"°*e'«" J2 ^me'^"% (98) 

m=-N 

with 



Am = ^l"^ F{z)e-'"-^^'dz. (99) 
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Obviously, one obtains only an approximation to the desired longitudinal pattern, because 
the trigonometric series (l98l) is necessarily truncated (A^ < A^max)- Its total number of 
terms, let us repeat, is fixed once the values of Q, L and ojq have been chosen. 
When p 7^ 0, the wavefield \E'(p, z,t) becomes 

N 

v|/(p,z,t) = e-^^o*e*«^ A^MkpmP)e'"-^"^\ (100) 

m=-N 

with 

= ^0- + • (101) 

The coefficients Am will yield the amplitudes and the relative phases of each Bessel beam 
in the superposition. 

Because we are adding together zero-order Bessel functions, we can expect a high field 
concentration around p = 0. Moreover, due to the known non-diffractive behaviour of the 
Bessel beams, we expect that the resulting wavefield will preserve its transverse pattern 
in the entire interval Q < z < L. 

The present methodology addresses itself to the longitudinal intensity pattern control. 
Obviously, we cannot get a total 3D control, due the fact that the field must obey the 
wave equation. However, we can use two ways to have some control over the transverse 
behaviour too. The first is through the parameter Q of eq. (p6!) . Actually, we have some 
freedom in the choice of this parameter, and FWs representing the same longitudinal 
intensity pattern can possess different values of Q. The important point is that, in su- 
perposition ( llOOp . using a smaller value of Q makes the Bessel beams to have a higher 
transverse concentration (because, on decreasing the value of Q, one increases the value 
of the Bessel beams transverse wave numbers), and this will reflect in the resulting fleld, 
which will present a narrower central transverse spot. The second way to control the 
transverse intensity pattern is using higher order Bessel beams, and we shall show this in 
Section 5.1.1. 

Now, let us present a few examples of our methodology. 
First example: 

Let us suppose that we want an optical wavefleld with Aq = 0.632 pm, i.e. with uoq = 
2.98 X 10^^ Hz, whose longitudinal pattern (along its z-axis) in the range < z < L is 
given by the function 
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[z - k){z - h 
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for k < z < L 



for k < z < k 



elsewhere 



(102) 



where k = L/5 — Azi2 and k = L/5 + Azi2 with Azu = L/50; while Z3 = L/2 — AZ34 and 
/4 = L/2 + A234 with A2;34 = L/10; and, at last, /s = 4L/5 — Azse and /e = 4L/5 + A2;56 
with Az56 = L/50. In other words, the desired longitudinal shape, in the range < z < L, 
is a parabolic function for /i < 2; < ^2, a unitary step function for /s < z < k, and again a 
parabola in the interval k z < k, being zero elsewhere (within the interval < z < L, 
as we said). In this example, let us put L = 0.2 m. 

We can then easily calculate the coefficients Am, which appear in superposition fllOOl) . 
by inserting eq. fll02p into eq. (!99l) . Let us choose, for instance, Q = 0.999 ci;o/c- This choice 
yields for m a maximum value iV^ax = 316, as one can infer from eq. fl97|) . Let us underline 
that one is not compelled to use just N = 316, but can adopt for any values smaller 
than it; more generally, any value smaller than that calculated via inequality fl97|l . Of 
course, when using the maximum value allowed for A^, one gets a better result. 



In the present case, let us adopt the value N = 30. In Fig.([3T]a) we compare the 
intensity of the desired longitudinal function F{z) with that of the Frozen Wave, "^{p = 
0, z, t), obtained from eq. (p8!l by adopting the mentioned value N = 30. 
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Figure 31: (a) Comparison between the intensity of the desired longitudinal function 
F{z) and that of our Frozen Wave (FW), \E'(p = 0,z,t), obtained from eq. fl98l) . The 
solid line represents the function F{z), and the dotted one our FW. (b) 3D-plot of the 
field-intensity of the FW chosen in this case by us. 
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One can verify that a good agreement between the desired longitudinal behaviour and 
our approximate FW is already got for = 30. The use of higher values for can only 
improve the approximation. Figure (l3T]b) shows the 3D-intensity of our FW, given by 



eq. fllOOl) . One can observe that this field possesses the desired longitudinal pattern, while 
being endowed with a good transverse localization. 

Second example (controlling the transverse shape too): 

We wish to take advantage of this example for addressing an important question: We can 
expect that, for a desired longitudinal pattern of the field intensity, by choosing smaller 
values of the parameter Q one will get FWs with narrower transverse width [for the 
same number of terms in the series entering eq. fllOOl) ]. because of the fact that the Bessel 
beams in eq. fllOOp will possess larger transverse wave numbers, and, consequently, higher 
transverse concentrations. We can verify this expectation by considering, for instance, 
inside the usual range < z < L, the longitudinal pattern represented by the function 



[z-h){z-h) 

F{z) = { ih-li? ^ (^03) 

elsewhere 



-4 ^ " for h<z<h 



with /i = L/2 — /S.Z and I2 = L/2 + Az. Such a function has a parabolic shape, with its 
peak centered at L/2 and with longitudinal width 2Az/\^. By adopting Aq = 0.632 /im 
(that is, ujq = 2.98 x 10^^ Hz), let us use superposition fllOOp with two different values of Q: 
We shall obtain two different FWs that, in spite of having the same longitudinal intensity 
pattern, possess different transverse localizations. Namely, let us consider L = 0.06 m and 
Az = L/lOO, and the two values Q = 0.999 c^o/c and Q = 0.995 tuo/c In both cases the 
coefficients Am will be the same, calculated from eq. (j99ll using this time the value = 45 
in superposition ( llOOp . The results are shown in Figs. (32a) and (32b). Both FWs have 
the same longitudinal intensity pattern, but the one with the smaller Q is endowed with 
a narrower transverse width. 

In this way, we can get some control on the transverse spot size through the pa- 
rameter Q. Actually, eq.f llOOl) . which defines our FW, is a superposition of zero-order 
Bessel beams, and, due to this fact, the resulting field is expected to possess a transverse 
localization around p = 0. Each Bessel beam in superposition f llOOp is associated with a 
central spot with transverse size, or width, Apm ~ 2.4/ kpm- On the basis of the expected 
convergence of the series fllOOp . we can estimate the width of the transverse spot of the 
resulting beam as being 

2.4 2.4 , , 

"1=0 Juof^ /c2 - 



which is the same value as that for the transverse spot of the Bessel beam with m = 
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Figure 32: (a) The Frozen Wave with Q = 0.999 c^o/c and = 45, approximately 
reproducing the chosen longitudinal pattern represented by eg. (11031) . (b) A different 
Frozen wave, now with Q = 0.995c(Jo/c (but still with = 45) forwarding the same 
longitudinal pattern. We can observe that in this case (with a lower value for Q) a higher 
transverse localization is obtained. 

in superposition (llOOp . Relation (I104p can be useful: Once we have chosen the desired 
longitudinal intensity pattern, we can even choose the size of the transverse spot, and use 
relation (11041) for evaluating the corresponding needed value of parameter Q. For a more 
detailed analysis concerning the spatial resolution and residual intensity of the Frozen 
Waves, we refer the reader to ref . |T 3Tj . 

The Frozen Waves, corresponding to zero group-velocity, are a particular case of the 
suhluminal Localized Waves. Actually, like in the Superluminal case, the (more orthodox, 
in a sense) subluminal LWs can be obtained by suitable superpositions of Bessel beams. 
They have been till now almost neglected, however, for the mathematical difficulties 
met in getting analytic expressions for them, difficulties associated with the fact that 
the superposition integral runs now over a finite interval. In Ref. |141j we have shown, 
by contrast, that one can indeed arrive at exact (analytic) solutions also in the case 
of general subluminal LWs, and both in the case of integration over the Bessel beams' 
angular frequency cj, and in the case of integration over their longitudinal wavenumber 
kz- We shall come back to this point in the following. 



10.1.1 Increasing the control on the transverse shape by using higher-order 
Bessel beams 



Here, we are going to argue that it is possible to increase even more our control on the 
transverse shape by using higher-order Bessel beams in our fundamental superposition 
(1 1001) . This new approach can be understood and accepted on the basis of simple and 
intuitive arguments, which are not presented here, but can be found in ref. |131j . A brief 
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description of that approach follows below. 

The basic idea is obtaining the desired longitudinal intensity pattern not along the 
axis p = 0, but on a cylindrical surface corresponding to p = p' > 0. To do that, we first 
proceed as before: Once we have chosen the desired longitudinal intensity pattern F{z), 
within the interval < ;2 < L, we calculate the coefficients as before, i.e., = 

(1/L) Jq F{z) exp{—i2TTmz/ L) dz, and fep^ = ^/cUq — [Q + 27rm/L)^. Afterwards, we 
just replace the zero-order Bessel beams Jo{kprnP), in superposition fllOOp . with higher- 
order Bessel beams, J^ikpmP), to get 

N 

v|/(p,z,t) = e-^'^°*e*«^ A^Jp{kp^p)e''^"'\ (105) 

m=-N 

From this result, and on the basis of intuitive argumentsLlSlJ, we can expect that the 
desired longitudinal intensity pattern, initially constructed for p = 0, will approximately 
shift to p = p', where p' represents the position of the first maximum of the Bessel function, 
i.e., the first positive root of the equation d Jp(fcp,m=o p)/dp) |p' = 0. 

By such a procedure, one can obtain very interesting stationary configurations of field 
intensity, as "donuts" , cylindrical surfaces, and much more. 

In the following example, we show how to obtain, e.g., a cylindrical surface of sta- 
tionary light. To get it, within the interval < z < L, let us first select the longitudinal 
intensity pattern given by eg. (11 03 p . with li = L/2 — Az and I2 = L/2 + Az, and with 
Az = L/300. Moreover, let us choose L = 0.05 m, Q = 0.998 cjq/c, and use = 150. 

Then, after calculating the coefficients Am by eg. (199!) . we use to superposition (11050 . 
choosing, in this case, p = 4. According to the previous discussion, one can expect the 
desired longitudinal intensity pattern to appear shifted to p' ~ 5.318/ kp „i=o = 8.47pm, 
where 5.318 is the value of kp^m=o P for which the Bessel function J4{kp^m=o p) assumes 

its maximum value, with kp^rn=o = \J^l ~ Q^- The figure [33] below shows the resulting 
intensity field. 

In Fig. (!33l a) the transverse section of the resulting beam for ^ = L/2 is shown. The 
transverse peak intensity is located at p = 7.75 pm, with a 8.5% difference w.r.t. the 
predicted value of 8.47 pm. Figure ([33]b) shows the orthogonal projection of the resulting 
field, which corresponds to nothing but a cylindrical surface of stationary light (or other 
fields). 

We can see that the desired longitudinal intensity pattern has been approximately 
obtained, shifted, as desired, from p = to p = 7.75 pm; and the resulting field resembles 
a cylindrical surface of stationary light with radius 7.75 pm and length 238 pm. Donut-like 
configurations of light (or sound) are also possible. 
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Figure 33: (a) Transverse section at z = L/2 of the considered higher-order FW. (b) 
Orthogonal projection of the three-dimensional intensity pattern of the same higher-order 
FW. 

10.2 Stationary wavefields with arbitrary longitudinal shape in 
absorbing media: An extension of the method. 

When propagating in a non-absorbing medium, the so-called nondiffracting waves main- 
tain their spatial shape for long distances. However, the situation is not the same when 
dealing with absorbing media. In such cases, both the ordinary and the nondiffracting 
beams (and pulses) will suffer the same effect: an exponential attenuation along the prop- 
agation axis. We shall present an extension [124] of the method given above, with the 
aim of showing that, through suitable superpositions of equal-frequency Bessel beams, it 
is possible to obtain even in absorbing media nondiffracting beams , whose longitudinal 
intensity pattern can assume any desired shape within a chosen interval < z < L of the 
propagation axis z. 

As a particular example, we are going to obtain new nondiffracting beams capable to 
resist the loss effects, maintaining amplitude and spot size of their central core for long 
distances. 

It is important to stress that the energy absorption by the medium continues to 
occur normally, but the new beams have an initial transverse field distribution, such to 
reconstruct (notwithstanding the presence of absorption) their central cores for distances 
considerably longer than the penetration depths of ordinary (nondiffracting or diffracting) 
beams. In this sense, the new method can be regarded as extending, for absorbing media, 
the self- reconstruction properties |142j that usual Localized Waves are known to possess 
in loss-less media. 

In the same way as for lossless media, we construct a Bessel beam with angular 
frequency u and axicon angle 6 in the absorbing materials by superposing plane waves, 
with the same angular frequency u, and whose wave vectors lie on the surface of a cone 
with vertex angle 6. The refractive index of the medium can be written as n{uj) = 
nji{uj) + ini{uj), quantity being the real part of the complex refraction index and rii the 
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imaginary one, responsible for the absorbtion effects. For a plane wave, the penetration 
depth 5 for the frequency uj is given by 5 = 1/a = c/2ujni, where a is the absorption 
coefficient. Therefore, a zero-order Bessel beam in dissipative media can be written as 
ip = Jo{kpp)exp{ii3z)exp{—iut) with P = n{Lu)u cos 9/ c = rinu cos9/c + iniu cos 9/c = 
Pr + iPh kp = n^uj sin 9/c + iniusm9/c = kp^ + ikpi, and so /c^ = n^a;^/c^ — Thus, 
it results ip = Jo((fcpR + ikpi)p)exp{iPiiz)exp{—iut)exp{—Piz), where /3r, /c^r are the 
real parts of the longitudinal and transverse wave numbers, and f3i, kpi are the imaginary 
ones, while the absorption coefficient of a Bessel beam with axicon angle 9 is given by 
ao = 2Pi = 2niujcos9/c, its penetration depth being 6g = 1/ag = c/2ujn\cos9. 

Due to the fact that kp is complex, the amplitude of the Bessel function Jo{kpp) 
starts decreasing from p = till the transverse distance p = l/2kpi, and afterwards it 
starts growing exponentially. This behaviour is not physically acceptable, but one must 
remember that it occurs only because of the fact that an ideal Bessel beam needs an 
infinite aperture to be generated. However, in any real situation, when a Bessel beam is 
generated by finite apertures, that exponential growth in the transverse direction, starting 
after p = l/2kpi, will not occur indefinitely, stopping at a given value of p. Let us moreover 
emphasize that, when generated by a finite aperture of radius R, the truncated Bessel 
beam [T3T] possesses a depth of field Z = R/ tan^, and can be approximately described 
by the solution given in the previous paragraph, ior p < R and z < Z. 

Experimentally, to guarantee that the mentioned exponential growth in the transverse 
direction does not even start, so as to meet only a decreasing transverse intensity, the 
radius R of the aperture used for generating the Bessel beam should be i? < l/2kpi. 
However, as noted by Durnin et al., the same aperture has to satisfy also the relation 
R > 27r/fcpR. From these two conditions, one can infer that, in an absorbing medium, 
a Bessel beam with just a decreasing transverse intensity can be generated only when 
the absorption coefficient is a < 2/A, i.e., if the penetration depth is 5 > A/2. The 
present method does refer to these cases, i.e., it is always possible to choose a suitable 
finite aperture size in such a way that the truncated versions of all solutions, including 
the general one given by eq.(lll), will not develop any unphysical behaviour. Let us now 
outline the method [141j . 

Let us consider an absorbing medium with the complex refraction index n{uj) = 
nji{uj) + ini{uj), and the following superposition of 2N + 1 Bessel beams with the same 
frequency u: 

N 

<f{p,z,t)= AraJA{K^m + ^kpiJp)e'^''-'e-'-'e-^^-\ (106) 

m=-N 

where the m are integer numbers, the Am are constant coefficients (yet unknown), quan- 
tities /3r„^ and fcpR^ (A^ and fcpi^) are the real (the imaginary) parts of the complex lon- 
gitudinal and transverse wave numbers of the m-th Bessel beam in superposition (11061) : 
the following relations being satisfied 
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^2.=-^^-/?^ (107) 



'1081 



where /^^ = /3r„ + /Cp„ = /CpR,„ + i/Cpi,„, with kp^^/kpi^ = n^/nj. 

Our goal is finding out the values of the longitudinal wave numbers and the 
coefficients Am in order to reproduce approximately, inside the interval < z < L (on 
the axis p = 0), a freely chosen longitudinal intensity pattern that we call |F(z)p. The 
problem for the particular case of lossless media[128,129], i.e., when nj = — >^ (3i^ = 0, 
was solved in the previous subsection. For those cases, it was shown that the choice 
(3 = Q + 2nm/L, with Am = F{z)exp{—i27imz/ L) / L dz can be used to provide 
approximately the desired longitudinal intensity pattern |F(2;)p in the interval < z < L, 
and, at the same time, to regulate the spot size of the resulting beam by means of the 
parameter Q. Such parameter, incidentally, can be also used to obtain large field-depths 
and moreover to inforce the linear polarization approximation to the electric field for the 
TE electromagnetic wave (see details in refs.[128,129]). 

However, when dealing with absorbing media, the procedure described in the last 
paragraph does not work, due to the presence of the functions exp{—Pi^z) in the super- 
position fll06p . because in this case that series does not became a Fourier series when 
p = 0. On attempting to overcome this limitation, let us write the real part of the 
longitudinal wave number, in superposition fll06p . as 



^ 27rm 

f^R^ = Q + -r- (109) 



with 



^ 27rm u! 
0<Q + —^<nn-. 110 
L c 

where inequality (11 101) guarantees forward propagation only, with no evanescent waves. 
In this way, the superposition fll06p can be written 

TV 

v^(p,z,t) = e--*e^Q^ Y: AmJoiikpn^+tkpiJp) e^'-^^e-^^-^ (111) 

m=-N 

where, by usin g eq. ([IOHD, we have f3i^ = {Q + 2'Km/ L)ni/n ^, a nd kp^ = kp^^ + ikpi^ 
is given by eq. fll07p . Obviously, the discrete superposition fillip could be written as a 
continuous one (i.e., as an integral over /3r„) by taking L — > oo, but we prefer the discrete 
sum due to the difficulty of obtaining closed-form solutions to the integral form. 

Now, let us examine the imaginary part of the longitudinal wave numbers. The 
minimum and maximum values among the (3i^ are (/3i)min = {Q — 2ttN / L)ni/n-^ and 
(A)max = {Q + '^T!-N/L)ni/nji, the central one being given by /?i = (A)m=o = Qni/riK- 
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With this in mind, let us evaluate the ratio A = [{j3i)max — (A)mm]//5i = 4:nN/LQ. Thus, 
when A << 1, there are no considerable differences among the various jSj^, because 
^ j3i holds for all m. In the same way, there are no considerable differences among 
the exponential attenuation factors, since exp{—Pi^z) ^ exp(— So, when p = the 
series in the r.h.s. of eq.f lllip can be approximately considered a truncated Fourier series 
multiplied by the function exp(— and, therefore, superposition fill II) can be used to 
reproduce approximately the desired longitudinal intensity pattern (on p = 0), 

within < z < L, when the coefficients Am are given by 



the presence of the factor exp^Pjz) in the integrand being necessary to compensate for 
the factors exp{—(3i^z) in superposition fillip . Since we are adding together zero-order 
Bessel functions, we can expect a good field concentration around p = 0. 

In short, we have shown in this Section how one can get, in an absorbing medium, 
a stationary wave-field with a good transverse concentration, and whose longitudinal 
intensity pattern (on p = 0) can approximately assume any desired shape \F{z)\'^ within 
the predetermined interval < z < L. The method followed above — let us resume — is a 
generalization of a previous one[128,129], and consists in the superposition in eq.f lllip of 
Bessel beams whose longitudinal wave numbers are individuated by the real and imaginary 
parts given in eqs.f llOQp and (11081) . respectively, while their complex transverse wave 
numbers are given by eq. fll07p . Finally, the coefficients of the superposition are given by 
eq.f lll2l) . The method is justified when AnN/LQ << 1: happily enough, this condition is 
satisfied in a great number of situations. 

Regarding the generation of these new beams, once we have an apparatus capable 
of generating a single Bessel beam, we can just use an array of such apparatuses to 
generate a sum of Bessel beams, with the appropriate longitudinal wave numbers and 
amplitudes/phases [as specified by our method], thus producing the desired final beam. 
For instance, we can use[128,129] a laser illuminating an array of concentric annular aper- 
tures (located at the focus of a convergent lens), with the appropriate radii and transfer 
functions in order to be able to yield both the required longitudinal wave numbers (once 
a value for Q has been chosen) and the coefficients An of the fundamental superposition 



10.2.1 Some Examples 

For generality's sake, let us consider a hypothetical medium in which a typical XeCl 
excimer laser (A = 308nm — = 6.12 x lO^^Hz) has a penetration depth of 5 cm; i.e. 
an absorption coefficient a = 20m~^, and therefore ui = 0.49 x 10^^. Besides this, let 
us suppose that the real part of the refraction index for this wavelength is = 1.5 and 
therefore n = tzr -|- ini = 1.5 -|- i 0.49 x 10~^. Note that the value of the real part of the 




(112) 



dniD. 
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refractive index is not so important for us, since we are dealing with monochromatic wave 
fields. 

A Bessel beam with uo = 6.12 x lO^^Hz and with an axicon angle 6 = 0.0141 rad 
(thus, with a transverse spot of radius 8.4 /xm), when generated by an aperture, say, of 
radius R = 3.5 mm, can travel in vacuum a distance equal to Z = R/ tan 6' = 25 cm 
while resisting the diffraction effects. However, in the material medium considered here, 
the penetration depth of this Bessel beam would be only Zp = 5 cm. Now, let us set forth 
two interesting applications of the present method. 

First Example: Almost Undistorted Beams in Absorbing Media. 

We can use the extended method to obtain, in the same medium and for the same wave- 
length, an almost undistorted beam capable of preserving its spot size and the intensity 
of its central core for a distance many times larger than the typical penetration depth of 
an ordinary beam (nondiffracting or not). To this purpose, let us suppose that, in the 
considered material medium, we want a beam (with uj = 6.12 x lO^^Hz) that maintains 
amplitude and spot size of its central core for a distance of 25 cm, i.e., a distance 5 times 
greater than the penetration depth of an ordinary beam with the same frequency. We 
can model this beam by choosing the desired longitudinal intensity pattern |F(2;)p (on 
p = 0), within < z < L, to be given by the function 



and by putting Z = 25 cm, with, for example, L = 33 cm. 

Now, one can use the Bessel beam superposition (11 111) to reproduce approximately 
the selected intensity pattern. Let us choose Q = 0.9999ci;/c for the /5r^ in eq. (11091) . 
and = 20 (notice that, according to inequality (1 11 01) . A^ could assume a maximum 
value of 158.) After having chosen the values of Q, L and A^, the values of the complex 
longitudinal and transverse wave numbers of the Bessel beams happen to be defined by 
relations ( ]109p . ( 1108^ and ( 1107^ . Eventually, we can have recourse to eq.( 1112p . and find 
out the coefficients Am of the fundamental superposition (lllip . that defines the resulting 
stationary wave-field. Let us just note that the condition 4ttN/LQ << 1 is perfectly 
satisfied in this case. 

In Fig.llMl a) we can see the 3D field- intensity of the resulting beam. One can observe 
that the field possesses a good transverse localization (with a spot size smaller than 
10 /im), and is capable of maintaining spot size and intensity of its central core till the 
desired distance (a better result could be reached by using a higher value of A^). 

It is interesting to note that at that distance (25 cm), an ordinary beam would have 
got its initial field-intensity attenuated 148 times. 

As we said above, the energy absorption by the medium continues to occur normally; 
the difference is that these new beams have an initial transverse field distribution sophis- 
ticated enough to be able to reconstruct (even in the presence of absorption) their central 




1 for 0<z<Z 
elsewhere. 
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Figure 34: (a) Three-dimensional field-intensity of the resulting beam, (b) The resulting 
beam, in an orthogonal projection and in logaritmic scale. 

cores, till a certain distance. For a better visualization of this field-intensity distribution 
and of the energy flux, we show in Fig. llMl b) the resulting beam, in an orthogonal pro- 
jection and in logarithmic scale. It appears clear that the energy comes from the lateral 
regions, in order to reconstruct the central core of the beam. On the plane z = 0, within 
the region p < R = 3.5 mm, there is an uncommon field intensity distribution, it being 
very dispersed instead of concentrated. This uncommon initial field intensity distribution 
is responsible for the construction of the central core of the resulting beam, and for its 
reconstruction all along the distance z = 25 cm. Due to absorption, the beam (total) 
energy, flowing through different z planes, is not constant; but the energy flowing in the 
beam spot area, and the beam spot size itself, are conserved till the distance (in this case) 
z = 25 cm. 

Second Example: Beams in absorbing media with a growing longitudinal field intensity. 

Let us consider again the previous hypothetical medium, in which an ordinary Bessel 
beam with 6 = 0.0141 rad and u = 6.12 x lO^^Hz would have a penetration depth of 
5 cm. We aim at constructing now a beam that, instead of possessing a constant core- 
intensity till the position z = 25 cm, presents on the contrary a (moderate) exponential 
growth of its intensity, till that distance {z = 25 cm). 

Let us assume we wish to get the following longitudinal intensity pattern |F(z)p, in 
the interval < z < L: 



with Z = 25cm and L = 33 cm. Using again Q = 0.9999a;/c, and = 20, we 




exp{z/Z) for < z < Z 
elsewhere , 



(114) 
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can proceed as in the first example, calculating the complex longitudinal and transverse 
wave numbers of the Bessel beams, and finally the coefficients Am of the fundamental 
superposition (11 111) . 

In Fig. (1351) we can see the 3D field-intensity of the resulting beam. One can observe 
that the field presents the desired longitudinal intensity pattern, with a good transverse 
localization (a spot size smaller than 10 /xm). 




Figure 35: Three-dimensional field-intensity of the resulting beam, in AN absorbing 
medium, with a growing longitudinal field intensity. 



Obviously, the amount of energy necessary to construct these new beams is greater 
than that necessary to generate an ordinary beam in a non-absorbing medium. And it is 
also clear that there is a limitation on the depth of field of these new beams. In the first 
example, for distances larger than 10 times the penetration depth of an ordinary beam, the 
field-intensity in the lateral regions would result to be higher than that at the core, and 
the field would loose the usual characteristics of a beam (transverse field concentration); 
not to speak of the greater energy demand. 



11 Subluminal Localized Waves (or Bullets) 

In this Section, abandoning for a while the subject of the so-called "Frozen Waves", 
we want to face the more general problem of obtaining, in a simple way, localized (non- 
diffractive) subluminal pulses as exact analytic solutions to the wave equations. } 141j These 
new ideal subluminal solutions, which propagate without distortion in any homogeneous 
linear media, will be here obtained for arbitrarily chosen frequencies and bandwidths. 



77 



avoiding in particular any recourse to the non-causal (backward moving) components 
that so frequently plague the previously known localized waves. The new solutions 
are suitable superpositions of — zeroth-order, in general — Bessel beams, which can be 
performed either by integrating w.r.t. the angular frequency a;, or by integrating w.r.t. 
the longitudinal wavenumber : Both methods are expounded in this review. The 
first one will appear to be powerful enough; we shall present the second method as well, 
however, since it allows dealing once more — from a different starting point — also with 
the limiting case of zero-speed solutions (and furnishes a new way, in terms of continuous 
spectra, for obtaining our Frozen Waves, so promising also from the point of view of 
applications). Some attention is moreover paid to the known role of Special Relativity, 
and to the fact that the localized waves are expected to be transformed one into the other 
by suitable Lorentz Transformations. At last, we briefly treat the case of non axially- 
symmetric solutions, in terms of higher order Bessel beams. The analogous pulses with 
intrinsic finite energy, or merely truncated, will be considered elsewhere. We keep fixing 
our attention especially on electromagnetism and optics: But let us repeat that results of 
the same kind are valid whenever an essential role is played by a wave-equation [like in 
acoustics, seismology, geophysics, elementary particle physics (as we verified even in the 
slightly different case of the Schroedinger equation), and also gravitation (for which we 
have recently got stimulating new results), and so on]. 



11.1 A foreword about the Sub luminal Localized Waves 

For self-consistency, let us repeat here the following considerations. For more than ten 
years, the so-called (non-diffracting) "Localized Waves" (LW), which are new solutions to 
the wave equations (scalar, vectorial, spinorial,...), are in fashion, both in theory and in 
experiment. In particular, rather well-known are the ones with luminal or Superluminal 
peak- velocity [143] : Like the so-called X-shaped waves (see |144l I145j and refs. therein; for 
a review, see, e.g., ref. |146j ). which are supersonic in Acoustics [117], and Superluminal in 
Electromagnetism (see [148] and refs. therein). 

As we know, since Bateman |149] and later on Courant & Hilbert[150], it was rec- 
ognized, e.g., that luminal LWs exist, which are solutions to the wave equations. More 
recently, some attention [9- 13] started to be paid to the (more "orthodox", as we said) 
suhluminal LWs too. Let us recall that all the LWs propagate without distortion — and 
in a self-reconstructive way[14-16] — in a homogeneous linear medium (apart from local 
variations): In the sense that their square magnitude keeps its shape during propagation, 
while local variations are shown only by its real, or imaginary, part. 

As in the Superluminal case, the subluminal LWs can be obtained by suitable super- 
positions of Bessel beams. |141] They have been till now almost neglected, as we know, 
for the mathematical difficulties met in getting analytic expressions for them, difficulties 
associated with the fact that the superposition integral runs over a finite interval. We 
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shall here re-address the question of such subluminal LWs, showing, by contrast, that one 
can indeed arrive at exact (analytic) solutions, both in the case of integration over the 
Bessel beams' angular frequency cj, and in the case of integration over their longitudinal 
wavenumber k^. 

As already claimed, the present work is devoted to the exact, analytic solutions: i.e., 
to ideal solutions. The corresponding pulses with finite energy, or truncated, will be 
presented elsewhere. 

Let us recall that, in the past, too much attention was not even paid to Brittingham's 
1983 paper [159J, wherein he had shown the possibility of obtaining pulse- type solutions to 
the Maxwell equations, which propagate in free space as a new kind of speed-c "solitons" . 
That lack of attention was partially due to the fact that Brittingham had been able neither 
to get correct finite-energy expressions for such "wavelets" , nor to make suggestions about 
their practical production. Two years later, however, Sezginer [160] was able to obtain 
quasi-nondiffracting luminal pulses endowed with a finite energy. Finite-energy pulses do 
no longer travel undistorted, as we do know, for an infinite distance, but can nevertheless 
propagate without deformation for a long field-depth, much larger than the one achieved 
by ordinary pulses like the gaussian ones: Cf., e.g., refs. [19-24] and refs. therein. 

Only after 1985 the general theory of LWs started to be extensively developed [25- 
31,2,6,3], both in the case of beams and in the case of pulses. For reviews, see for instance 
the refs. |146[ I174[ 11651 11631 1166] and citations therein. For the propagation of LWs in 
bounded regions (like wave-guides), see refs. [33-36] and refs therein. For the focusing of 
LWs, see the Second Part of this review [as well as refs. [17911180] and quotations therein]. 
As to the construction of general LWs propagating in dispersive media, see refs. [39-47]; 
and, for lossy media, cf. ref. |158] and refs. therein. Al last, for finite-energy, or truncated, 
solutions see refs. [48-50,24,3,34], and work in progress. 

By now, the LWs have been experimentally produced [147l 11931 1194] . and are being 
applied in fields ranging from ultrasound scanningj l95l [T96| 1153] to optics (for the produc- 
tion, e.g., of new type of tweezers [197] ). All those works have demonstrated by now that 
nondiffr acting pulses can travel with an arbitrary peak-speed v, that is, with < f < oo; 
while Brittingham and Sezginer had confined themselves to the luminal case {v = c) only. 

As we were remarking, the Superluminal and luminal LWs have been, and are being, 
intensively studied; whilst the subluminal ones have been neglected: Almost all the few 
papers dealing with them had recourse till now to the paraxial [T98j approximation [199] . 
or to numerical simulations [154] . due to the above mentioned mathematical difficulty in 
obtaining exact analytic expressions for subluminal pulses. Indeed, only one analytic 
solution was known[9-ll, 28, 57,58], biased by the physically unconvenient facts that its 
frequency spectrum is very large, that it doesn't even possess a well-defined central fre- 
quency, and, even more, that backward-travelling [168 1 IT66j components (ordinarily called 
"non-causal", since they should be entering the antenna or generator) were needed for 
constructing it. Aim of the next Sections is showing, on the contrary, that subluminal 
localized exact solutions can be constructed with any spectra, in any frequency bands 
and for any bandwidths; and without employing [145[ 1165] any backward-travelling com- 
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ponents. 



12 A first method for constructing physically 
acceptable, subluminal Localized Pulses 

Axially symmetric solutions to the scalar wave equation are known to be superpositions of 
zero-order Bessel beams over the angular frequency uj and the longitudinal wavenumber 
kz'- That is, in cylindrical co-ordinates, 



= dcuy_^^^dA;,5(a;,A;,)Jo|py— -A;2je*^^^e-*-*, (115) 



where /c^ = tu^/c^ — k1 is the transverse wavenumber. Quantity fc^ has to be positive 
since evanescent waves cannot come into the play. 

The condition characterizing a nondiffr acting wave is the existence [166], |20T] of a linear 
relation between longitudinal wavenumber and frequency uo for all the Bessel beams 
entering superposition (113); that is to say, the chosen spectrum has to entail p^Sl 1163] 
for each Bessel beam a linear relationship of the typeF*! 



oj = vkz + b (116) 

with 6 > 0. Requirement (116) can be regarded also as a specific space-time coupling, 
implied by the chosen spectrum S. Equation (116) has to be obeyed by the spectra of any 
one of the three possible types (subluminal, luminal or Superluminal) of nondiffr acting 
pulses. Let us mention that with the choice in eq.(116) the pulse re-gains its initial 
shape after the space-interval Azi = 2'Kv/h. But the more general case can be also 
considered |145l 1189] when h assumes any values bm = rrib (with m an integer), and the 
periodicity space- interval becomes Azm = Azi/m . We are referring ourselves, of course, 
to the real (or imaginary) part of the pulse, since its modulus is known to be endowed 
with rigid motion. 

In the subluminal case, of interest here, the only exact solution known till recent time, 
represented by eq.(124) below, was the one found by Mackinnon |151j . Indeed, by taking 
into explicit account that the transverse wavenumber kp of each Bessel beam entering 
eq.(115) has to be real, it can be easily shown (as first noticed by Salo et al. for the 
analogous acoustic solutions [151]) that in the subluminal case b cannot vanish, but must 
be larger than zero: b > 0. Then, on using conditions (116) and 6 > 0, the subluminal 
localized pulses can be expressed as integrals over the frequency only: 

**More generally, as shown in ref. [145j . the chosen spectrum has to call into the play, in the plane 
u, kz, if not exactly the line (114), at least a region in the proximity of a straight-line of that a type. It 
is interesting that in the latter case one obtains solutions endowed with finite energy, but possessing a 
finite "depth of field" , that is to say, nondiffracting only till a certain finite distance. 
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\E'(p, = exp[—ib-] / du S{uj) Jo{pkp) exp[iLJ- 



where now 



with 



kp = -^2bu-b^-{l-vyc^)iu^ 



(117) 



:il8) 



C = z-vt 



and with 



1 + v/c 



:il9) 



(120) 



1 — v/c 

As anticipated, the Bessel beam superposition in the subluminal case resuhs to be an in- 
tegration over a finite interval of u, which does clearly shows that the backward-travelling 
(non-causal) components correspond to the interval < u < b. It could be noticed 
that eq.(117) does not represent the most general exact solution, which on the contrary 
is a sum \189\ of such solutions for the various possible values of b mentioned above: That 
is, for the values b^ = mb with spatial periodicity A^;^ = I\zxjm. But we can confine 
ourselves to solution (117) without any real loss of generality, since the actual problem 
is evaluating in analytic form the integral entering eq.(117). For any mathematical and 
physical details, see ref. [189j . 

Now, if one adopts the change of variable 



equation (117) becomes |154j 



1 — jc? 



:i2i] 



^{p,z,t) 



exp [—■i-z\ exp 

V 




ds S{s) J( 




(122) 



In the following we shall adhere — as it is an old habit of ours — to some symbols standard 
in Special Relativity (since the whole topic of subluminal, luminal and Superluminal LWs 
is strictly connected [HSl EM EDS] with the principles and structure of SR [cf.[203l EQlj 
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and refs. therein], as we shall mention also in the concluding remarks which fllow below); 
namely: 



As already said, eq.(122) has till now yielded one analytic solution for S{s) = 
constant, only (for instance, S{s) = 1); which means nothing but S{uj) = constant: 
in this case one gets indeed the Mackinnon solution \151\ \170\ \153\ 1192] 



c 



i- P'y'^ T] 
c 



xsinc J^72(p2 + 7'C') , (124) 



c2 



which however, for its above-mentioned drawbacks, is endowed with little physical and 
practical interest. In eq.(124) the sine function has the ordinary definition 



sine X = (sin x)/x 

and 

7] = z-Vt, with V= — , (125) 



where V and v are related by the de Broglie relation. [Notice that \l/ in eq.(124), and in 
the following ones, is eventually a function (besides of p) of z, t via ( and rj, both functions 
of z and t ]. 

However, we can construct by a very simple method new subluminal pulses corre- 
sponding to whatever spectrum, and devoid of backward- moving (i.e., "entering") com- 
ponents, just by taking advantage of the fact that in our equation (122) the integration 
interval is finite: that it, by transforming it in a good, instead of a harm. Let us first ob- 
serve that eq.(122) doesn't admit only the already known analytic solution corresponding 
to S{s) = constant, and more in general to ^(ti;) = constant, but it will also yield an 
exact, analytic solution for any exponential spectra of the type 

S{uj) = exp [-^] , (126) 

with n any integer number, which means for any spectra of the type S{s) = 
exp[m7r//3] exp[m7rs], as can be easily seen by checking the product of the various 
exponentials entering the integrand. In eq.(126) we have set 
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Q = a;+ — a;_ 



The solution writes in this more general case: 



*(p,C,?7) = 26/37' exp 



i- rj 
c 



X exp [in—\ smc 



62 (h ^ 

_ ^2 p2 _^ - 7^ ( + nTT 



(127) 



Let us exphcitly notice that also in eq.(127) quantity r] is defined as in Eqs.(125) above, 
where V and v obey the de Broglie relation vV = c', the subluminal quantity v being 
the velocity of the pulse envelope, and V playing the role (in the envelope's interior) of a 
Superluminal phase velocity. 

The next step, as anticipated, consists just in taking advantage of the finiteness of 
the integration limits for expanding any arbitrary spectra S{u!) in a Fourier series in the 
interval cu^ < cu < cu^ : 



°° 271 

S{uj) = ^ ^„exp[+m— a;] 



where (we went back, now, from the s to the uu variable): 

1 2n 
= ^ y dujS{uj) exp[-m— cj] 



(128) 



(129) 



quantity fl being defined as above. 

Then, on remembering the special solution (127), we can infer from expansion (126) 
that, for any arbitrary spectral function S{lj), one can work out a rather general axially- 
symmetric analytic solution for the subluminal case: 



^^(^,(,77) = 26/3 7' exp 



i- (5^'^ r] 
c 



X An exp [irJ^] sine 

n=—oo P 



i 



62 (h ^ 



(130) 



in which the coefficients An are still given by eq.(129). Let us repeat that our solution is 
expressed in terms of the particular equation (127), which is a "Mackinnon-type" solution. 

The present approach presents many advantages. We can easily choose spectra local- 
ized within the prefixed frequency interval (optical waves, microwaves, etc.) and endowed 
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with the desired bandwidth. Moreover, as we have seen, spectra can now be chosen 
such that they have zero value in the region U-^ < uj < b, which is responsible for the 
backward-travelling components of the subluminal pulse. 

Let us stress that, even when the adopted spectrum S{uj) does not possess a known 
Fourier series (so that the coefficients An cannot be exactly evaluated via eq.(129)), one 
can calculate approximately such coefficients without meeting any problem, since our 
general sohitions (130) will still be exact solutions. 

Let us set forth some examples. 



12.1 Some Examples 

In general, optical pulses generated in the lab possess a spectrum centered at some 
frequency value, ujq, called the carrier frequency. The pulses can be, for instance, ultra- 
short, when Acu/cuo > 1; or quasi- monochromatic, when Acu/cuo << 1, where Acu is the 
spectrum bandwidth. 

These kinds of spectra can be mathematically represented by a gaussian function, or 
functions with similar behaviour. 

First Two Examples: 

Let us first consider a gaussian spectrum 



whose values are negligible outside the frequency interval uj_ < oo < ujj^ over which the 
Bcsscl beams superposition in eq.(117) is made, it being uj_ = 6/(1 + (3) and uj^ = 
6/(1 — P). Of course, relation (116) has still to be satisfied, with 6 > 0, for getting an 
ideal subluminal localized solution. Notice that, with spectrum (131), the bandwidth 
(actually, the FWHM) results to be Aa; = 2/a. Let us emphasize that, once v and b 
have been fixed, the values of a and uq can afterwards be selected in order to kill the 
backward-travelling components, that correspond, as we know, to lu < b . 

The Fourier expansion in eq.(128), which yields, with the above spectral function 
(131), the coefficients 



constitutes an excellent representation of the gaussian spectrum (131) in the interval 
< u < (provided that, as we requested, our gaussian spectrum does get negligible 
values outside the frequency interval a;_ < a; < a;+). 




(131) 



A ~ 




(132) 
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In other words, on choosing a pulse velocity v < c and a value for the parameter 6, 
a subluminal pulse with the above frequency spectrum (131) can be written as eq.(129), 
with the coefficients An given by eq.(132): the evaluation of such coefficients An being 
rather simple. Let us repeat that, even if the values of the An are obtained via a (rather 
good, by the way) approximation, we based ourselves on the exact solution eq.(130). 

One can, for instance, obtain exact solutions representing subluminal pulses for optical 
frequencies. Let us get the subluminal pulse with velocity v = 0.99 c, carrier angular 
frequency ujq = 2.4 x 10^^ Hz (that is, Aq = 0.785 /xm) and bandwidth (FWHM) Au = 
a;o/20 = 1.2 x 10^^ Hz, which is an optical pulse of 24 fs (which is the FWHM of the 
pulse intensity). For a complete pulse characterization, one has to choose the value of 
the frequency b: let it be 6 = 3 x 10^'^ Hz; as a consequence one has u!_ = 1.507 x 10^^ Hz 
and uj^ = 3 X 10^^ Hz. [This is exactly a case in which the considered pulse is not 
plagued by the presence of backward-travelling components, since the chosen spectrum 
possesses totally negligible values for u < b]. The construction of the pulse does already 
result satisfactory when considering about 51 terms (—25 < n < 25) in the series entering 
eq.(130). 

Figures [36] show our pulse, evaluated by summing the mentioned fifty-one terms. 
Namely: Fig. (36a) depicts the orthogonal projection of the pulse intensity; Fig. (36b) 
shows the three-dimensional intensity pattern of the real part of the pulse, which reveals 
the carrier wave oscillations. 




(a) P] 



Figure 36: (a) The intensity orthogonal projection for the pulse corresponding to 
eqs.(131,132) in the case of an optical frequency (see the text); (b) The three-dimensional 
intensity pattern of the real part of the same pulse, which reveals the carrier wave oscil- 
lations. 

Let us stress that the ball-like shap^ for the field intensity should be typically asso- 
ciated with all the subluminal LWs, while the typical Superluminal ones are known to be 

ttit can be noted that each term of the series in eq.(130) corresponds to an elhpsoid or, more specifically, 
to a spheroid, for each velocity v. 
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X-shaped [144^ I2U2] . as predicted since long by special relativity in its "non-restricted" 
version: See refs. |203l 12041 11481 1146j and refs therein. 

A second spectrum ^(ti;) would be, for instance, the "inverted parabola" one, centered 
at the frequency Uq: that is. 



(133) 

otherwise , 



where Au, the distance between the two zeros of the parabola, can be regarded as the 
spectrum bandwidth. One can expand S{u), given in eq.(133), in the Fourier series (128), 
for uj^ < uj < UJ+ , with coefficients that — even if straightforwardly calculable — results 
to be complicated, so that we skip reporting them here explicitly. Let us here only mention 
that spectrum (133) may be easily used to get, for instance, an ultrashort (femtasecond) 
optical non- diffracting pulse, with satisfactory results even when considering very few 
terms in expansion (128). 

Third Example: 

As a third interesting example, let us consider the very simple case when — within the 
integration hmits U-, uj+ — the complex exponential spectrum (124) is replaced by the 
real function (still linear in uj) 

S{uj) = ^ -r exp[a{uj -UJ+] , (134) 

1 — exp [—0(^7+ — 

with a a positive number [for a = one goes back to the Mackinnon case]. Spectrum 
(134) is exponentially concentrated in the proximity of u;+, where it reaches its maximum 
value; and becomes more and more concentrated (on the left of a;_|_, of course) as the 
arbitrarily chosen value of a increases, its frequency bandwidth being Aa; = 1/a. Let us 
recall that, on their turn, quantities 0;+ and (x;_ depend on the pulse velocity v and on 
the arbitrary parameter b. 

By performing the integration as in the case of spectrum (126), instead of solution 
(127) in the present case one eventually gets the solution 



X exp 



2abP'y^ exp [067^] exp [—au^ 



c 



1 — exp [—a{uj^ — u;_)] 

b 

c 



smc 



(135) 
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After Mckinnon's, this eq.(135) appears to be the simplest closed-form solution, since 
both of them do not need any recourse to series expansions. In a sense, our solution (135) 
may be regarded as the subluminal analogue of the (Superluminal) X-wave solution; a 
difference being that the standard X-shaped solution has a spectrum starting with 0, 
where it assumes its maximum value, while in the present case the spectrum starts at u;_ 
and gets increasing afterwards, till uj+. More important is to observe that the gaussian 
spectrum has a priori two advantages w.r.t. eq.(134): It may be more easily centered 
around any value ujq of a;, and, when increasing its concentration in the surroundings of 
tuo, the spot transverse width does not increase indefinitely, but tends to the spot-width 
of a Bessel beam with uj = ujq and = [ujq ~ b)/V, at variance with what happens for 
spectrum (134). Anyway, solution (135) is noticeable, since it is really the simplest one. 

Figure [37] shows the intensity of the real part of the subluminal pulse corresponding 
to this spectrum, with v = 0.99 c, with 6 = 3 x 10^^ Hz (which result in uj_ = 1.5 x 10^^ Hz 
and 0;+ = 3 X 10^^ Hz), and with Au/u^ = 1/100 (i.e., a = 100). This is an optical pulse 
of 0.2 ps. 



( Re '■]' f 




Figure 37: The intensity of the real part of the subluminal pulse corresponding to spectrum 
(134), with V = 0.99 c, with 6 = 3 x 10^^ Hz (which result in U- = 1.5 x 10" Hz and 
= 3 X 10^^ Hz), and with Acj/cj+ = 1/100 (i.e., a = 100). 
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13 A second method for constructing 
subluminal Localized Pulses 



The previous method appears to be very efficient for finding out analytic subluminal LWs, 
but it looses its validity in the limiting case v ^ 0, since for v = it is U- = a;+ and the 
integral in eq.(117) degenerates, furnishing a null value. By contrast, we are interested 
also in the v = case, since it corresponds, as we already know, to some of the most 
interesting, and potentially useful, LWs: Namely, to the stationary solutions to the wave 
equations endowed with a static envelope, and that we have called "Frozen Waves" . 

Before going on, let us recall that the theory of Frozen Waves was initially developed 
in refs. |191l 1205] . by having recourse to discrete superpositions in order to bypass the 
need of numerical simulations. In the case of continuous superpositions, some numerical 
simulations were performed in refs. |206j . However, the method presented in this Section 
does allow finding out analytic exact solutions (without any further need of numerical 
simulations) even for Frozen Waves consisting in continuous superpositions. Actually, we 
are going to see that the present method works whatever is the chosen field-intensity 
shape, also in regions with size of the order of the wavelength. 

It is possible to get such results by starting again from eq.(115), with constraint (114), 
but going on — this time — to integrals over k^, instead of over uj. It is enough to write 
relation (116) in the form 

k, = -(u-b) (116') 

V 

for expressing the exact solutions (115) as 

*(p,2;,t) = exp[-ibt] / dk^ S{k^) Jo{pkp) exp[iCk^] , (136) 



with 



/e,mm = — 



kzinax 



c l + (3 
b 1 



(135) 



and with 



cl-p 



kp' = -^ + 2-(3k, + -, (138) 
7^ c c 

where quantity ( is still defined according to eq.(119), always with v < c. 

One can show that the unique exact solution previously known [151] may be rewritten 
in form (136) with S{kz) = constant. Then, on following the same procedure exploited in 
our first method (previous Section), one can find out new exact solutions corresponding 
to 
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S{h) = exp[— ^] , (139) 

where 



K = /c^max — /c^min , 

by performing the change of variable [analogous, in its finality, to the one in eq.(121)] 

h = -^^s + l3). (140) 
c 

At the end, the exact subluminal solution corresponding to spectrum (139) results to 

be 

b b 
*(p,C,r;) = 2-7' exp[i-pj^r]] 



X exp [iriTT/?] sine 



%i^p'+(-l'( + nA , (141) 



We can again observe that any spectra S{kz) can be expanded, in the interval /c^min < 
max, in the Fourier series: 

S{k,) = An exp [+tn—k,] , (142) 

n=—oo 

with coefficients given now by 

^ 1? ^kzS{kz) ex^[-in—kz] (143) 

i\ ^fczinin -A. 

quantity having been defined before. 

At the end of the whole procedure, the general exact solution representing a sublu- 
minal LW, for any spectra S{kz), can be eventually written 

b b 
c c 



00 

X 'Y An exp [iriTT/S] sine 

n=— 00 



^'7^+ (-'yH + riTr] , (144) 



\ c 



whose coefficients are expressed in eq.(143), and where quantity 77 is defined as above, in 
eq.(125). 

Interesting examples could be easily worked out, as we did at the end of the previous 
Section. 
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14 Stationary solutions with zero-speed envelopes 



Here, we shall refer ourselves to the (second) method, expounded in the previous Section. 
Our solution (144), for the case of envelopes at rest, that is, in the case v = [which 
imphes ( = z], becomes 



b °° 
\l'(p, z,t) =2 - exp[—ibt] ^ sine 



\ 




(145) 



with coefficients An given by eq.(143) with f = 0, so that its integration limits simplify 
into —b/c and b/c, respectively. Thus, one gets 



= ^ y_,//^^ ^^^^^ [-myM . (143') 



Equation (145) is a new exact solution, corresponding to stationary beams with a 
static intensity envelope. Let us observe, however, that even in this case one has energy 
propagation, as it can be easily verified from the power fiux Sg = — V\E'7^ d'^n/dt (scalar 
case) or from the Poynting vector = {E A H) (vectorial case: the condition being that 

be a single component, Az, of the vector potential A). |148j We have indicated by 
the real part of ^ . For f = 0, eq.(116) becomes 

uj = b = uq , 

so that the particular subluminal waves endowed with null velocity are actually monochro- 
matic beams. 

It may be stressed that the present (second) method does yield exact solutions, with- 
out any need of the paraxial approximation, which, on the contrary, is so often used when 
looking for expressions representing beams, like the gaussian ones. Let us recall that, 
when having recourse to the paraxial approximation, the obtained beam expressions are 
valid only when the envelope sizes (e.g., the beam spot) vary in space much more slowly 
than the beam wavelength. For instance, the usual expression for a gaussian beamfT98] 
holds only when the beam spot Ap is much larger than Aq = uJq/{2'kc) = 6/(27rc): so that 
those beams cannot be very much localized. By contrast, our method overcomes such 
problems, since we have seen that it yields exact expressions for (well localized) beams 
with sizes of the order of their wavelength. Notice, moreover, that the already-known 
exact solutions — for instance, the Bessel beams — are nothing but particular cases of our 
solution (145). 

An example: On choosing (with < g_ < g+ < 1) the spectral double-step function 
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uJo{q+ - q- 







for q_ujQ/c < kz < q+Uo/c 
elsewhere , 



(146) 



the coefficients of eq.(145) become 



2nnuJo{q+ - g_) 



(147) 



The double-step spectrum (146) corresponds, with regard to the longitudinal wave 
number, to the mean value = uJo{q+ + q-)/'^c and to the width Ak^ = uJo{q^ — q^)/c. 
From such relations, it follows that Akz/k^ = 2(g+ — q_)/{q+ + qj). 

For values of g_ and g+ that do not satisfy the inequahty Ak-^jk^ « 1, the resulting 
solution will be a non-paraxial beam. 

Figures 38 show the exact solution corresponding to uiq = 1.88 x 10^^ Hz (i.e., Aq = 
1 /im), and to g_ = 0.3, and to g+ = 0.9: It results to be a beam with a spot-size 
diameter of 0.6 /im, and, moreover, with a rather good longitudinal localization. In the 
case of Eqs.(144, 145), about 21 terms (—10 < n < 10) in the sum entering eq.(143) are 
quite enough for a good evaluation of the series. The beam considered in this example 
is highly non-paraxial (with Ak^/kz = 1 ), and therefore couldn't have been obtained by 
ordinary gaussian beam solutions (which are valid in the paraxial regime only)!^]. 




Figure 38: (a) Orthogonal projection of the three-dimensional intensity pattern of the 
beam (a null-speed subluminal wave) corresponding to spectrum (146); (b) 3D plot of 
the field intensity. The beam considered in this example is highly non-paraxial. 



ttWe are considering here only scalar wave fields. In the case of non-paraxial optical beams, the vector 
character of the field has to be taken into account 
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14.1 A new approach to the "Frozen Waves" 

Let us now emphasize that a noticeable property of the present method is that it 
allows a spatial modeling even of monochromatic fields (that correspond to envelopes 
at rest] so that, in the electromagnetic cases, one can speak, e.g., of the modeling of 
"light-fields at rest"). 

Let us recall that such a modelling — rather interesting, especially for 
applications [197] — was already performed in refs. |191l 12051 1158j . and has been exploited 
at the beginning of our Third Part [cf. eqs.(99, 100)], in terms of discrete superpositions 
of Bessel beams. And the stationary fields with static envelopes have been called "Frozen 
Waves" (FW) by us. 

But the method presented in the last Sections allows us to make use of continuous 
superpositions, in order to get a predetermined longitudinal (on-axis) intensity pattern, 
inside a desired space interval < z < L. In fact, the continuous superposition, analogous 
to eq.(lOO), now writes 



which is nothing but the previous eq.(136) with v = (and therefore ( = z). In other 
words, eq.(148) does just represent a null-speed subluminal wave. To be clearer, let us 
recall once more that the FWs were expressed in the past as discrete superposition because 
it was not known at that time how to treat analytically a continuous superposition like 
(148). We are now able, however, to extend the previous approach to FWs to the case 
of integrals: without numerical simulations, but in terms once more of analytic solutions. 

Indeed, the exact solution of eq.(148) is given by eq.(145), with coefficients (143'); 
and one can choose the spectral function S{kz) in such a way that assumes the on-axis 
pre-chosen static intensity pattern |F(2;)p. Namely, the equation to be satisfied by S{kz), 
to such an aim, comes out by associating eq.(148) with the requirement |\E'(p = 0, z, t)\^ = 
\F{z)\'^, which entails the integral relation 



Equation (149) would be trivially solvable in the case of an integration between — oo and 
+00, since it would merely be a Fourier transformation; but obviously this is not the 
case, because its integration limits are finite. Actually, there are functions F{z) for which 
eq.(149) is not solvable at all, in the sense that no spectra S{kz) exist obeying the last 
equation. For instance, if we consider the Fourier expansion 




(148) 




(149) 



F{z) = / dkz S{kz) & 



J — oo 



when ^'(/cz) does assume non-negligible values outside the interval —ujq/c < k^ < ujq/c, 
then in eq.(149) no S{kz) can forward that particular F{ result. 
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However, way-outs can be devised, such that one can nevertheless find out a function 
S{kz) that approximately (but satisfactorily) complies with eq.(149). 

The first way-out consists in writing S{kz) in the form 

Sih) = ^ f: F (^) e-^^"-^^/^ , (150) 

where, as before, K = 2ujq/c. Then, one can easily verify eq.(150) to guarantee that the 
integral in eq.(149) yields the values of the desired F{z) at the discrete points z = 2x171 /K. 
Indeed, the Fourier expansion (150) is already of the same type as eq.(142), so that in this 
case the coefficients of our solution (145), appearing in eq.(143'), do simply become 

1 , 2n7r, 

An = F{-^) . (151) 

This is a powerful way for obtaining a desired longitudinal (on-axis) intensity pattern, 
especially for tiny spatial regions, because it is not necessary to solve any integral to find 
out the coefficients An, which by contrast are given directly by eq.(151). 

Figures [39] depict some interesting applications of this method. A few desired longi- 
tudinal intensity patterns |F(^)p have been chosen, and the corresponding Frozen Waves 
calculated by using eq.(145) with the coefficients An given in eq.(151). The desired pat- 
terns are enforced to exist within very small spatial intervals only, in order to show the 
capability of the method to model p^T] the field intensity shape also under such strict 
requirements. 

In the four examples below, we considered a wavelength A = 0.6 //m, which corre- 
sponds to Wo = = 3.14 X 10^^ Hz. 

The first longitudinal (on-axis) pattern considered by us is that given by 



Fiz) 



^a{z-Z) {]< z< Z 



elsewhere 



i.e., a pattern with an exponential increase, starting from z = Q untill z = Z. The chosen 
values of a and Z are Z = 10 fim and a = 3/Z. The intensity of the corresponding Frozen 
Wave is shown in Fig. (39a). 

The second longitudinal pattern (on-axis) taken into consideration is the gaussian 
one, given by 



e-^d )' for -Z <z< Z 



Fiz) 



elsewhere , 
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with q = 2 and Z = l.Q /iin. The intensity of the corresponding Frozen Wave is shown in 
Fig. (39b). 

In the third example, the desired longitudinal pattern is supposed to be a super- 
gaussian: 



«(f)^™ for -Z<z<Z 



F(z) 



elsewhere , 



where m controls the edge sharpness. Here we have chosen q = 2, m = 4 and Z = 2 fim. 
The intensity of the Frozen Wave obtained in this case is shown in Fig. (39c). 

Finally, in the fourth example, let us choose the longitudinal pattern as being the 
zero-order Bessel function 



F{z) 



Jo{qz) for — Z < z < Z 



elsewhere 



with g = 1.6 X 10^ m~^ and Z = 15 fim. The intensity of the corresponding Frozen Wave 
is shown in Fig.(39d). 

Let us observe that any static envelopes of this type can be easily transformed into 
propagating pulses by the mere application of Lorentz transformations. 

Another way-out exists for evaluating S{kz), based on the assumption that 

S{h) ~ S{h) , (152) 

which consitutes a good approximation whenever S{kz) assumes negligible values outside 
the interval [—wq /c, ujo/c]. In such a case, one can have recourse to the method associated 
with eq.(142) and expand S{kz) itself in a Fourier series, getting eventually the relevant 
coefficients An by eq.(143). Let us recall that it is still K = fc^max — fcj-min = 2uoq/c. 

It is worthwhile to call attention to the circumstance that, when constructing FWs 
in terms of a sum of discrete superpositions of Bessel beams (as it has been done by us 
in the Second part of this work, and in refs. |191l 12051 11581 1197j ). it was easy to obtain 
extended envelopes like, e.g., "cigars": Where "easy" means by using only a few few terms 
of the sum. By contrast, when we construct FWs — following this Section — as continuous 
superpositions, then it is easy to get highly localized (concentrated) envelopes. Let us 
explicitly mention, moreover, that the method presented in this Section furnishes FWs 
that are no longer periodic along the 2;-axis (a situation that, with our old method [191l 
12051 1158j . was obtainable only when the periodicity interval tended to infinity). 
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Figure 39: Frozen Waves with the on-axis longitudinal field pattern chosen as: (a) Expo- 
nential; (b) Gaussian; (c) Super-gaussian; (d) Zero order Bessel function 



15 The role of Special Relativity, and of Lorentz 
Transformations 

Strict connections exist between, on one hand, the principles and structure of Special 
Relativity and, on the other hand, the whole subject of subluminal, luminal, Superluminal 
Localized Waves, in the sense that it is expected since long time that a priori they are 
transformable one into the other via suitable Lorentz transformations (cf . refs. |203[ 12041 
120 7j . besides work of ours in progress).^ 

Let us first confine ourselves to the cases faced in this Third Part. Our subluminal 
localized pulses, that may be called "wave bullets", behave as particles: Indeed, our 

tt'^ Let us call attention to a paper by Saari et al.[2D8j, noticed by us only recently, wherein the 
relativistic connections among the LWs were also investigated in terms of suitable LTs. We are actually 
glad in quoting here such a reference because it appears inspired by a philosophy, which — going back in 
part to papers like refs. j^03l I204j — has been constantly shared by our old, and recent, papers. Let us 
mention also a further interesting article noticed by us recently, by Besieris et al. [209| . wherein ordinary 
LTs were already used, correctly, in the context of subluminal LWs [whilst the Superluminal LTs used in 
that paper seem, however, to be partially defective]. 
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subluminal pulses [as well as the luminal and Superluminal (X-shaped) ones, that have 
been amply investigated in the past literature] do exist as solutions of any wave equations, 
ranging from electromagnetism and acoustics or geophysics, to elementary particle physics 
(and even, as we discovered recently, to gravitation physics). From the kinematical point 
of view, the velocity composition relativistic law holds also for them. The same is true, 
more in general, for any localized waves (pulses or beams). 

Let us start for simplicity by considering, in an initial reference-frame O, just a {u- 
order) Bessel beam: 

^(p, 0, z, t) = Mpkp) e'"^ e''''^ e"*"* . (153) 

In a second reference-frame O', moving with respect to (w.r.t.) O with speed u — along the 
positive z-axis and in the positive direction, for simplicity's sake — , it will be observed |208j 
the new Bessel beam 

^{p',(f)\z',t') = Mp'k'p,) e""^' e'''^'^' e-'"'*' , (154) 
obtained by applying the appropriate Lorentz transformation (a Lorentz "boost") with 

k' p' = kp] k'z' = 'y{kz — uuj / oj' = 7(0; — ukz) ; (155) 
this can be easily seen, e.g., by putting 

p = p'; z = -f{z' + ut'); t = -i{t' + uz' /c^) (156) 

directly into eq.(154). 

Let us now pass to subluminal pulses. We can investigate the action of a Lorentz 
transformation (LT), by expressing them either via the first method (Section 12) or via 
the second one (Section 13). Let us consider for instance, in the frame O, a f-speed 
(subluminal) pulse, given by eq.(155) of our Section 15. When we go on to a second 
observer O' moving with the same speed v w.r.t. frame O, and, still for the sake of 
simplicity, passing through the origin O of the initial frame at time t = 0, the new 
observer O' will see the pulse [208] 

^{p,z',t') = e-'*'"'° n du S{uj) Jo{p'k'p>) e''''''^' , (157) 

Ju)- 

with 

k' z' = ^~^uj /v — ^h/v; Lu' = 'yb; k' p/ = lu'q/c^ — k' z/'^ , (158) 

as one gets from the Lorentz transformation in eq.(155), or in eq.(156), with u = v [and 
7 given by Eqs.(123)]. Notice that k'z' is a function of u, as expressed by the first one of 
the three relations in the previous Eqs.(158); and that here uj' is a constant. 
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If we explicitly insert into eq.(157) the relation uj = 'y{vk'z' + 7&), which is nothing 
but a re- writing of the first one of Eqs.(158), then eq.(157) becomes |208] 



■LUo/c 



.iz'k' ^1 



(159) 



where S is expressed in terms of the previous function S{ijj)^ entering eq.(157), as follows: 



Equation (159) describes monochromatic beams with axial symmetry (and does coincide 
also with what derived within our second method, in Section 13, when posing = 0). 

The remarkable conclusion is that a subluminal pulse, given by our eq.(117), which 
appears as a f-speed pulse in a frame O, will appear ^208j in another frame O' (travelling 
w.r.t. observer O with the same speed v in the same direction z) just as the monochromatic 
beam in eq.(159) endowed with angular frequency oj'q = 76, whatever be the pulse spectral 
function in the initial frame O: even if the kind of monochromatic beam, one arrives to, 
does of course depend^ on the chosen S{uj). The vice- versa is also true, in general. 

Let us set forth explicitly an observation that up to now has been noticed only in 
ref. |141j . Namely, let us mention that, when starting not from eq.(117) but from the 
most general solutions which — as we have already seen — are sums of solutions (117) 
over the various values bm of b, then a Lorentz transformation will lead us to a sum of 
monochromatic beams: actually, of harmonics (rather than to a single monochromatic 
beam). In particular, if one wants to obtain a sum of harmonic beams, one has to apply 
a LT to more general subluminal pulses. 

Let us add that also the various Superluminal localized pulses get transformed |208j 
one into the other by the mere application of ordinary LTs; while it may be expected that 
the subluminal and the Superluminal LWs are to be linked (apart from some known tech- 
nical difficulties, that require a particular caution) by the Superluminal Lorentz "transfor- 
mations" expounded long ago, e.g., in refs. |204[ I210[ 12071 1203] and refs. therein. Let us 
recall once more that, in the years 1980-82, special relativity, in its non-restricted version, 
predicted that, while the simplest subluminal object is obviously a sphere (or, in the limit, 
a space point), the simplest Superluminal object is on the contrary an X-shaped pulse (or, 
in the limit, a double cone): cf. FigJTT| taken^ from refs. |203| I204j . The circumstance 

tt+ One gets in particular a Bessel-type beam when 5' is a Dirac's delta-function: S{uj) = 6{uj — ujq). 
Moreover, let us notice that, on applying a LT to a Bessel beam, one obtains another Bessel beam, with 
a different axicon-angle. 

W++ One should pay attention that, as we were saying, the topic of Superluminal LTs is a delicate [204|, 
l210U207ll203| one, at the extent that the majority of the recent attempts to re-address this question and 
its applications seem to be defective (sometimes they do not even keep the necessary covariance of the 
wave equation itself). 

tt* Let us recall, more specifically, that Fig. 11 depicts the following. Let us start from an object that be 
intrinsically spherical, i.e., that is a sphere in its rest-frame (Panel (a)). Then, after a generic subluminal 
LT along x, i.e., under a subluminal x-boost, it is predicted by Special Relativity (SR) to appear as 




(160) 
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that the locahzed solutions to the wave equations follow the same behaviour is rather 
interesting, and is expected to be useful — in the case, e.g., of elementary particles and 
quantum physics — for a deeper comprehension of de Broglie's and Schroedinger's wave 
mechanics. With regard to the fact that the simplest subluminal LWs, solutions to the 
wave equation, are "ball-like" , let us present in FigsHUl in ordinary 3D space, the general 
shape of the Mackinnon's solutions, as expressed by eq.(124) for v << c: In such figures 
we graphically depict the field iso-intensity surfaces, which result to be (as expected) just 
spherical in the considered case. 

We have also seen, among the others, that, even if our first method (Section 12) 
cannot directly yield zero-speed envelopes, such envelopes "at rest", eq.(145), can be 
however obtained by applying a f-speed LT to eq.(130). In this way, one starts from 
many frequencies [eq.(130)] and ends up with one frequency only [eq.(145)], since b gets 
transformed into the frequency of the monochromatic beam. 



16 Non-axially symmetric solutions: The case of 
higher-order Bessel beams 

Let us stress that till now we paid attention to exact solutions representing axially- 
symmetric (subluminal) pulses only: that is to say, to pulses obtained by suitable super- 
positions of zero-order Bessel beams. 

It is however interesting to look also for analytic solutions representing non-axially 
symmetric subluminal pulses, which can be constructed in terms of superpositions of u- 
order Bessel beams, with u a positive integer (z/ > 0). This can be attempted both in 
the case of Sect. 12 (first method), and in the case of Sect. 13 (second method). 

For brevity's sake, let us take only the first method (Sect. 12) into consideration. One 
is immediately confronted with the difficulty that no exact solutions are known for the 

ellipsoidal due to Lorentz contraction (Panel (b)). After a Super luminal a:-boost [2041 1210|, I207j (namely, 
when this object moves with Superluminal speed V), it is predicted by SR — in its non-restricted 
version (ER) — to appear [203j as in Panel (d), i.e., as occupying the cylindrically symmetric region 
bounded by a two-sheeted rotation hyperboloid and an indefinite double cone. The whole structure, 
according to ER, is expected to move rigidly and, of course, with the speed V, the cotangent square of 
the cone semi-angle being (V/c)'^ — 1. Panel (c) refers to the limiting case when the boost-speed tends to 
c, either from the left or from the right (for simplicity, a space axis is skipped). It is remarkable that the 
shape of the localized (subluminal and Superluminal) pulses, solutions to the wave equations, appears 
to follow the same behaviour; this can have a role for a better comprehension even of de Broglie and 
Schroedinger wave-mechanics. See also Fig. 40. 
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Figure 40: In the previous Figure we have seen how SR, in its non-restricted version (ER), 
predicted [2031 [201] that, while the simplest snhlummal object is obviously a sphere (or, in 
the hmit, a space point), the simplest Superluminal object is on the contrary an X-shaped 
pulse (or, in the limit, a double cone). The circumstance that the Localized Solutions to 
the wave equations do follow the same pattern is rather interesting, and is expected to 
be useful — in the case, e.g., of elementary particles and quantum physics — for a deeper 
comprehension of de Broglie's and Schroedinger's wave mechanics. With regard to the 
fact that the simplest subluminal LWs, solutions to the wave equations, are "ball-like" , let 
us depict by these figures, in the ordinary 3D space, the general shape of the Mackinnon's 
solutions as expressed by eq.(124), numerically evaluated for v << c. In figures (a) and 
(b) we graphically represent the field iso-intensity surfaces, which in the considered case 
result to be (as expected) just spherical. 

integral in eq.(122) when Jo(.) is replaced with Ju{-)- 

One can overcome this difficulty by following a simple method, which allows obtaining 
"higher-order" subluminal waves in terms of the axially-symmetric ones. Indeed, it is 
well-known that, if \l/(a;, 2;, t) is an exact solution to the ordinary wave equation, then 
d'^/dx and d'^/dy are also exact solutions.* By contrast, when working in cylindrical 
co-ordinates, if \l/(p, 0, 2, t) is a solution to the wave equation, quantities d'^jdp and 
d'^/d(f) are not solutions, in general. Nevertheless, it is not difficult at all to reach the 
noticeable conclusion that, once '${p,(f), z,t) is a solution, then also 

*(p,,,M)^e-*(^ + --) (161) 

tt* Let us mention that even and will be exact solutions. 
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is an exact solution! For instance, for an axially-symmetric solution 

of the type = Jo{kpp) exp[ikz] exp[—iujt], equation (11611) yields \E' = 

—kp Ji{kpp) exp[z0] exp[ifc2] exp[— zcjt], which is actually one more analytic solution. 

In other words, it is enough to start for simplicity from a zero-order Bessel beam, 
and to apply eg. (11611) . successively, v times, in order to get as a new solution ^ = 
{—kpY Ju{kpp) exp[iz/0] exp[i/c2] exp[— zcut], which is a i/-order Bessel beam. 

In such a way, when applying v times eq. (ll6ip to the (axially-symmetric) subluminal 
solution \E'(p, 2;,t) in Eqs. (130,129,128) [obtained from eq.(117) with spectral function 
S{uj)], we get the subluminal non-axially symmetric pulses \E'v(p, 0, ^, t) as new analytic 
solutions, consisting as expected in superpositions of z^-order Bessel beams: 

*„(p, 0, z,t) = / duo S\uj) Mkpp) e'"^ e''^' e'^' , (162) 

where kp{uj) is given by eq.(118), and quantities S'^uj) = {—kp{u)y S{uj) are the spectra 
of the new pulses. If S{u!) is centered at a certain carrier frequency (it is a gaussian 
spectrum, for instance), then S'{uj) too will approximately result to be of the same type. 

Now, if we wish the new solution \E',^(p, 0, z, t) to possess a pre-defined spectrum 
S'{u!) = F^uj), we can first take eq.(117) and put S{uj) = F{ui) / {—kp{ijj)Y in its solution 
(130), and afterwards apply to it, v times, the operator U = exp[i0] [d/dp+ {i / p)d / dcj))] : 
As a result, we will obtain the desired pulse, \I'j,(p, 0, 2;, t), endowed with 5"(ci;) = F{ijj). 

An example: 

On starting from the subluminal axially-symmetric pulse \E'(p, 2;, t), given by eq.(130) 
with the gaussian spectrum (131), we can get the subluminal, non-axially symmetric, 
exact solution \I'i(p, 0, z, t) by simply calculating 

vl/i(p,0,^,t) = 1^ (163) 

which actually yields the "first-order" pulse \E'i(p, 0, 2;, t), which can be more compactly 
written in the form: 



^i(p, 0,^,0 = 2-t;72 exp 
c 



c 



An exp [in-] (164) 



with 



where 



fe2 

,2 



V'in(p, 0, VX) = -J 7V Z-^ [Z COS Z - Sin Z] e"^ , (165) 



c 



\ 




^7V+ ( ^7'C + ^7r) . (166) 
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This exact solution, let us repeat, corresponds to superposition fll62p . with S'{uj) = 
kp{(jj)S{ijj), quantity ^(ci;) being given by eq. fll3ip . It is represented in Figure IHl The 
pulse intensity has a "donut-like" shape. 




Figure 41: Orthogonal projection of the field intensity corresponding to the higher order 
subluminal pulse represented by the exact solution eq.(163), quantity \l/ being given by 
eq.(128) with the gaussian spectrum (131). The pulse intensity happens to have this time 
a "donut"-like shape. 



16.1 A few concluding remarks about the Third Part 

In this Third Part we started developing, by suitable superpositions of equal-frequency 
Bessel beams, a first theoretical and experimental methodology to obtain localized sta- 
tionary wave fields, with high transverse localization, whose longitudinal intensity pattern 
can approximately assume any desired shape within a chosen interval < 2; < L of the 
propagation axis z. Their intensity envelope remains static, i.e., with velocity f = 0; so 
that we named "Frozen Waves" (FW) these new solutions to the wave equations (and, 
in particular, to the Maxwell equations). Inside the envelope of a FW only the carrier 
wave does propagate: And the longitudinal shape, within the interval < 2; < L, can be 
chosen in such a way that no nonnegligible field exists outside the pre-determined region 
(consisting, e.g., in one or more high intensity peaks). Such solutions are noticeable also 
for the different and interesting applications they can have, especially in electromagnetism 
and acoustics, such as optical tweezers, atom guides, optical or acoustic bistouries, various 
important medical apparatus (mainly for destroying cancerous cells), etc. 

Afterwards, we have addressed the more general subject of the subluminal Localized 
Waves, and shown that — like in the well-known Superluminal case |143j — the subluminal 
solutions can be obtained by superposing Bessel beams [T4T] . Such solutions have been 
scarcely considered in the past, for the reason that the superposition integral has to 
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run in this case over a finite interval (which makes mathematically difficult to work out 
analytic expressions for them). We have shown, however, how it is possible to obtain, 
in a simple way, non- diffracting subluminal pulses as exact analytic solutions to the wave 
equations: For arbitrarily chosen frequencies and bandwidths, and avoiding any recourse 
to the backward-travelling components. 

Indeed, till recent times only one closed-form subluminal LW solution, ipcf, to the 
wave equations was known [TST] : obtained by choosing in the relevant integration a con- 
stant weight-function S{uj)] whilst all other solutions had been previously got only by 
numerical simulations. On the contrary, a subluminal LW can be obtained in closed form 
by adopting, for instance, any spectra S{uj) that be expansions in terms of ipci- In fact, 
the initial disadvantage, of having to deal with a limited bandwidth, may be turned into 
an advantage, since in the case of "truncated" integrals the spectrum ^(ti;) can be ex- 
panded in a Fourier series. More in general, it has been shown how can one arrive at 
exact solutions both by integration over the Bessel beams' angular frequency u, and by 
integration over their longitudinal wavenumber k^. Both methods have been expounded 
above. The ffist one appears to be comprehensive enough; we have studied the second 
method as well, however, since it furnishes a new way, in terms of continuous spectra, 
for tackling also the limiting case of zero-speed solutions (i.e., for obtaining the Frozen 
Waves) . 

We have briefly treated the case, moreover, of non axially- symmetric solutions, that 
is, of higher order Bessel beams. 

At last, some attention has been paid to the role of Special Relativity, and to the fact 
that the localized waves are expected to be transformable one into the other by suitable 
Lorentz Transformations. Moreover, our results seem to show that in the subluminal case 
the simplest LW solutions are (for v << c) "ball"-like, as expected since long |203j on the 
mere basis of special relativity |204j . [Indeed, let us recall once more that already in the 
years 1980-82 it had been predicted that, if the simplest subluminal object is a sphere (or, 
in the limit, a space point), then the simplest Superluminal object is an X-shaped pulse 
(or, in the limit, a double-cone); and viceversa: Cf. Figs. 11. It is rather interesting that 
the same pattern appears to be followed by the localized solutions of the wave equations]. 
For the subluminal case, see, e.g.. Figs. 40. 

The subluminal localized pulses, endowed with a finite energy, or merely truncated, 
will be presented elsewhere. 
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